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volumes  are  the  results  of  the  model  and  these  may  be  directly  input  to  the 
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optional  output  of  the  model. 
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r,KCTION  I 


INTRODUCTION 


1. 1  Background 

The  objective  of  this  project  was  to  develop  an  air  base  motor 
vehicle  model  that  would  simulate  an  air  base  traffic  system  by  using 
available  transportation  and  land  use  data.  The  model  produces  a  graphi¬ 
cal  representation  of  the  air  base  motor  vehicle  operation  and  also 
produces  a  file  of  traffic  flows  for  input  to  the  Air  Quality  Assessment 
Model  (AQAM) ,  Thus,  AQAIl  will  be  made  sensitive  to  air  base  land  use  and 
transportation  activities,  and  will  be  able  to  predict  motor  vehicle- 
related  air  pollution  effectively. 

The  Base  Automotive  Transportation  Simulation  (BATS)  Model  is 
useful  for  assessing  the  effectiveness  of  various  planning  strategies 
such  as  strategies  to  reduce  and  control  air  pollution,  traffic  engi¬ 
neering  projects  to  reduce  congestion  and  improve  flow,  energy  conserva¬ 
tion  strategies,  and  land  use  strategies  represented  in  Air  Force  com¬ 
prehensive  plans. 

The  model  was  tested  and  its  performance  evaluated  at  three  typical 
air  bases  (Tinker,  Williams,  and  Davis-Monthan) .  Data  for  inputs  to  the 
model  were  collected  by  the  project  team  in  two  phases:  The  first  phase 
included  a  survey  and  extraction  of  existing  data  from  the  TAB  A-1  and 
Military  Traffic  Management/Transportation  Engineering  Agency  data.  The 
second  phase  involved  field  surveys  at  the  three  air  bases.  Information 
was  collected  on  base  operations  necessary  to  operate  the  model  (e.g., 
work  schedules  and  the  peak  operational  periods  of  such  facilities  as 
gas  stations,  hospitals,  and  base  exchanges). 

The  data  collection  task  is  documented  in  a  technical  report  that 
provides  specific  guidance  to  air  base  personnel  who  use  the  model  (see 
ESL-TR-79-16,  Volume  I,  Data  Collection  and  Reduction).  The  data  collec¬ 
tion  document  demonstrates  step  by  step  how  to  extract  required  data 
from  the  existing  sources  and  how  to  collect  data  by  using  field  survey 
techniques,  if  required.  The  document  also  provides  average  or  default 
values  (which  can  be  used  as  model  inputs)  that  are  based  on  tlie  three 
air  base  studies  undertaken  by  the  project  team. 


1. 2  Review  of  Regulations 

On  February  25,  1974  (Federal  Register,  Vol.  39,  p.  7270),  the 
Environmental  Protection  Agency  administrator  promulgated  a  regulation 
(40  CFR  52.22  b)  for  the  review  of  indirect  sources  for  Inclusion  in  the 
Clean  Air  Act  implementation  plans.  Subsequent  changes  and  clarifications 


1 


were  issued  (Federal  Register.  Vol.  39.  p.  25292,  July  9,  1974)  and  the 
regulations  were  scheduled  for  implementation  on  January  1,  1975.  As  a 
result  of  congressional  action  in  December  1974,  the  EPA  was  prohibited 
from  using  FY  1975  funds  for  the  enforcement  of  indirect-source  (and  so- 
called  parking  management)  regulations. 

According  to  the  regulations,  indirect  sources  include,  but  are  not 
limited  to:  highways  and  roads;  parking  facilities;  retail,  commercial, 
industrial,  educational,  recreational,  amusement,  sports,  and  entertain¬ 
ment  facilities;  office  and  government  buildings;  apartment,  condominium, 
and  housing  projects;  and  airports.  For  those  facilities  that  may  be 
aggregated  as  parking  or  trip  attraction  facilities,  the  regulations 
require  that  an  air  quality  Impact  assessment  be  made  if  the  new  or  incre¬ 
mental  parking  capacity  exceeds  certain  limits.  Within  Standard  Metro¬ 
politan  Statistical  Areas  (SMSA) ,  this  includes  new  facilities  with 
parking  capacity  of  1000  cars  or  more  or  an  increase  of  an  existing 
facility  by  500  cars  or  more.  In  the  case  of  highways,  the  regulations 
apply  if  the  expected  average  dally  traffic  (ADT)  meets  or  exceeds  20,000 
within  ten  years  for  a  new  roadwav,  or  If  there  is  an  increase  of  10,000 
or  more  on  a  modified  roadway.  Outside  the  SMSAs ,  all  the  limits  are 
exactly  double. 

To  satisfy  all  technical  requirements  for  the  assessment  of  the 
anticipated  air  quality  impact  of  all  roadway  and  parking-related  projects, 
the  Air  Force  has  developed  the  comprehensive  AQAM  and  BATS  models.  These 
models  provide  a  detailed  air  quality  impact  analysis  by  first  performing  a 
dynamic  analysis  of  traffic  flow  within  and  adjacent  to  the  base  and  then 
determining  emissions  and  atmospheric  dispersion  of  pollutants. 


In  addition  to  providing  a  comprehensive  analysis  tool  for  assess¬ 
ment  of  the  local  impact  of  roadways  and  parking  facilities  (in  terms  of 
ambient  concentrations  of  carbon  monoxide),  BATS  and  AQAM  can  be  used  for 
a  variety  of  similar  analytical  requirements.  First,  they  can  be  used  to 
meet  the  requirements  for  preparing  environmental  impact  statements  under 
the  National  Environmental  Policy  Act  of  1970  (Public  Law  91-190)  and 
similar  state  laws.  BATS  and  AQAM  would,  of  course,  also  be  directly 
applicable  to  satisfying  the  local  impact-assessment  requirements  of  the 
proposed  (and  now  delayed)  parking  management  regulations  (Federal 
Register,  Vol.  39,  p.  30440,  August  22,  1974). 


I . 3  Computer  Requirements 

BATS  is  written  in  FORTRAN  IV.  It  has  been  debugged  and  tested  on 
SRI's  CDC  6400  computer  and  the  CDC  6600  at  Eglin  AFB.  It  has  been  run 
successfully  under  the  K.R0N0S  and  SCOPE,  operating  systems  using  60,000 
words  of  core  sLf)rage.  The  model  should  be  readily  transferable  to  any 
computer  system  witli  FORTRAN  IV  capabilitv*  Becausi-  of  thi-  large  core 
:.l  I'l  .U'.i’  I  e<|ii  I  I  emeiil  s  ,  howevii,  the  ptogi'am  woiiltl  need  mod  1  I  I  ('<  1 1  I  o  u  lielore 
I  I  coulil  he  lull  on  nii'dinm  oi  cm.ill  '.i-.ile  coiiipnli'l  '.vsleiiis.  Jnicli  mod  I  I  I 
cations  could  he  nnderlaken  vi  (lioni  miii.li  dlMii  ultv  so  tliat  the  IIAI'S 
program  could  lie  run  using  32  kilobites  of  computer  storage. 
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1. 4  Document  OrRanlzation 

This  document  describes  the  BATS  model  and  discusses  each  subroutine 
or  program  in  the  model.  Section  2  provides  the  logical  functioning  of 
the  model,  the  basic  equations  used  in  the  model,  and  the  reasons  behind 
the  unique  features  of  the  model.  Therefore,  the  second  section  is  use¬ 
ful  to  the  person  who  wants  to  know  the  general  features  of  the  model. 

Section  3  is  useful  to  the  person  who  wants  to  understand  and 
perhaps  change  an  operational  result  of  the  model.  This  section 
discusses  the  purpose  of  each  subroutine  or  program,  its  Inputs,  and 
its  outputs.  A  listing  of  each  routine  is  provided  in  Appendix  A. 
Appendix  B  describes  the  utility  routines.  Appendix  C  presents  a  sample 
run  of  BATS  for  Williams  AFB,  and  Appendices  D  and  E  (Volume  3)  contain 
the  sample  runs  of  Tinker  AFB  and  Davis- Monthan  AFB,  respectively. 
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SKC'I  ION  II 


BATS— (JKNKKAL  DESCRIPTION 


2.1  Introduction 


BATS  combines  the  features  of  transportation  planning  and  traffic 
flow  analysis  by  predicting  volumes,  speeds  and  delays  on  the  network  of 
streets  and  zones  which  comprise  an  air  base.  The  model  synthesizes  the 
capabilities  of  the  SRI  Network  Analysis  Program  (SNAP)  (Haney  and 
Thompson,  1971)  and  Indirect  Source  Model  for  Pollution  (ISMAP)  (Sandys 
et  al. ,  1975). 

Figure  1  presents  a  flow  diagram  of  the  BATS  program  at  its  most 
general  level.  The  model  is  very  straightforward,  consisting  of  INPUT, 
INITIALIZATION,  TRANSPORTATION  PLANNING,  TRAFFIC  FLOW,  and  OUTPUT  func¬ 
tions.  Each  iteration  causes  a  loop  through  the  entire  program.  Because 
this  loop  includes  the  input  function,  the  user  inputs  new  data  for  each 
hour  to  be  run.  In  this  way,  the  number  of  trips  and  other  network 
characteristics  may  be  changed  for  any  time  period  simulated. 

All  tests  and  loops  shown  in  the  flowchart  are  required  to  perform 
multiple  iterations  through  the  model.  Multiple  iterations  allow  the 
user  flexibility  in  the  accuracy,  amount  of  computer  inputs,  and  computer 
time  associated  with  each  computer  run  of  the  model.  For  example,  the 
user  may  select  1  hour  as  the  length  of  time  during  which  constant  demand 
and  flow  will  provide  desired  simulation  accuracy.  He  may  also  choose 
to  simulate  12  hours  of  the  day  as  representing  most  traffic  flow  on  the 
base.  Then  the  program  would  input  12  sets  of  data  and  perform  12  itera¬ 
tions  to  predict  traffic  during  12  hours.  Alternately,  if  a  user  requires 
less  accuracy,  he  might  choose  to  simulate  12  hours  in  one  iteration; 
therefore,  the  time  period  simulated  would  be  12  hours  and  only  one  set 
of  input  data  would  be  read. 

The  flow  diagram  also  shows  a  special  feature  that  allows  for  15- 
minute  iterations  without  looping  through  the  input  function.  Typically, 
congestion  will  not  occur  during  an  entire  hour  on  a  street  network; 
rather,  it  is  confined  to  a  shorter,  perhaps  15-minute,  time  period. 
Therefore,  the  user  should  use  the  15-minute  Iteration  when  he  wishes  to 
predict  possible  congestion  occurring  on  the  base.  Wlien  les:;  accuracy 
is  required,  a  1-,  2-,  12-,  or  24-hour  time  period  may  suffice  to  pre¬ 
dict  travel  on  the  network.  The  output  of  the  model  includes  printed 
reports,  plots  of  the  traffic  network,  and  a  traffic  data  file  for  input 
to  the  AQAM  program.  The  output  is  described  in  the  BATS  User  Guide 
ESI,-TR-79- 1 6  Volume  I:  Data  (Collection  and  Reduction  and  therefore  is 
Mol  Inchrrled  in  this  document. 

/  I  p  I  p  1 1 1  , 

Dala  -PIP-  ipppppi  Ippi  p.ipIp  I  1  pppi  pp  i  I  od  ni  1  I  e  P  .p  I  I  ppip  Ih.ll  hp  p.pil  ol  ;i 
'  "ipppppl  p-r  tpppp.  iipp'  1/  1  ypi's  ipI  jppppil  ,irp-  idetPl  i  I  i  P'd  by  G.prd  Types  (( 
t  Ip  rpppp;',lp  I  h .  G.pril',  ppI  I  yp<-  (I  .iiiil  I  .are  re.pd  I  (pr  e.ich  i  t  e  r.p  t  i  pup  . 
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INPUT  -  READS  ALL  DATA  INPUT 
REQUIRED  FOR  ONE  ITERATION 
OF  THE  MODEL,  (USUALLY  DATA 
INPUTS  FOR  1  HOUR) 


INITIALIZE  -  DETERMINES  THE  SIZE  OF 
EACH  PARKING  ZONE,  PLOTS  A  LINK  MAP 
IF  CALLED  FOR,  DETERMINES  TRAVEL 
TIME  OVER  EACH  LINK  IN  THE  NETWORK 


TRANSPORTATION  PLANNING  FUNCTIONS  - 
GENERATES  TRIPS  BETWEEN  ZONES, 
DISTRIBUTES  TRIPS  BETWEEN  ZONE  PAIRS, 
DETERMINES  THE  MODE  OF 
TRANSPORTATION,  AND  ASSIGNS  VEHICLES 
TO  THE  TRAFFIC  NETWORK 


TRAFFIC  FLOW  ANALYSIS  -  DETERMINES 
INTERSECTION  QUEUEING  AND  DELAY, 
DETERMINES  LINK  DELAY,  DETERMINES 
PARKING  LOT  RUNNING  TIME  AND  FLOW 


OUTPUT  -  PRINTS  RESULTS  OF  EACH 
MODEL  FUNCTION,  PLOTS  VOLUME-FLOW 
MAPS  OF  THE  BASE,  WRITES  AQAM  PARAMETERS 
OF  SPEED  AND  VOLUME  FOR  INPUT  TO  AQAM 


They  describe  the  heading,  the  options,  and  the  number  of  Card  Types  2 
through  16  to  be  read  for  an  iteration.  Card  Types  2,  3,  4,  5,  and  10, 
which  define  the  network,  need  to  be  read  only  on  the  first  Iteration 
of  the  program.  The  network  definition  includes  link,  intersection,  zone, 
demographic  variable,  and  truck  route  specifications.  Card  Types  6,  7, 

8,  9,  and  11  contain  information  that  may  vary  from  one  time  period  to 
another,  such  as  trip  generation  coefficients,  gate  counts,  through 
traffic  counts,  load  factors,  and  work-shift  counts.  Card  Types  12,  13, 
14,  15,  and  16  are  special  purpose  cards,  which  reset  default  values  of 
plot  scale  factors,  calibration  coefficients,  default  trip  generation 
arrays,  demographic  variable  names,  and  daily-to-yearly  traffic  volume 
conversion  factors. 

The  input  routine  is  designed  to  minimize  input  cards  by  permanently 
storing  values  read  on  a  previous  iteration,  and  using  these  if  no  other 
values  are  read. 


2.3  I  n  i  t  Lai  i ^  t_l  on 

The  second  function  shown  on  the  flowchart  stores  the  initial  values 
of  parameters  that  are  required  by  the  simulation.  Three  parameters  are 
initiall zed : 

•  The  x-y  coordinates  of  each  link  are  read  by  INPUT  and  oriented 
by  the  direction  of  traffic  flows. 

•  The  area  of  each  zone  is  used  to  compute  the  average  length  of 
each  trip  within  the  zone  needed  to  find  a  parking  spot.  This 
length  is  the  longest  side  of  the  zone,  plus  the  area  divided  by 
the  longest  link. 

•  The  time  to  travel  each  link  is  determined  based  on  the  speed 
limit  and  length  of  each  link.  On  later  iterations,  the  travel 
time  is  a  function  of  the  link  and  intersection  capacity  and  flow. 

2 . 4  Transportation  Planning 

The  third  function  shown  on  the  flowchart  predicts  the  network 
vehicle  counts  using  the  transportation  planning  methodology.  A  great 
many  different  "transportation  planning"  methodologies  could  have  been 
used  to  accomplish  this  function.  A  study  for  small  urban  areas  (Grecco 
et  al.,  1976),  undertaken  as  part  of  the  National  Cooperative  Highway 
kesearch  Program,  is  a  useful  reference  defining  several  of  the  approaches. 

An  air  base  is  quite  similar  to  a  small  urban  area  in  that  it  has 
housing,  office,  industry,  recreation,  shopping,  and  service  facilities. 
However,  an  air  base  dii fers  from  many  urban  areas  because  a  great  deal 
of  Information  has  been  collected  about  the  various  facilities.  The 
TAB  A-1  Environmental  Narrative  found  at  most  air  bases  Illustrates  the 
availability  of  data  that  can  be  used  in  the  TRANSPORTATION  PLANNING 
p  rocess . 
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2.4.1  General  Me 


'I'Ik-  |)ii  r|)C);;c  of  Liu-  t  ranspor  Lot  i  on  planning  mctliO(loloj<y  is  to  <-st  i  fti;it< 
lor  selected  Lime  peri<jds,  tfie  nuinber  of  velii(  les  on  ea(  fi  link  of  a 
network.  The  methodology  proceeds  through  five  sequential  steps  as  shown 
In  the  flowchart.  Figure  2.  The  first  step,  trip  generation  (TKIPGKN) , 
is  the  most  complex  step  in  the  methodology;  it  allows  the  user  many 
optional  ways  of  generating  trips  from  or  to  each  zone.  Trip  generation 
predicts  how  many  persons  wish  to  make  a  trip  away  from  (productions)  or 
to  (attractions)  a  zone.  The  origin  or  the  destination  of  a  trip  is  pre¬ 
dicted,  but  origin  and  destination  pairs  must  be  predicted  by  a  subse¬ 
quent  function. 

The  trip  distribution  function  (GRAVO)  determines  the  number  of 
person  trips  going  between  each  origin-destination  (0-D)  pair.  In  this 
function,  trips  that  go  from  an  exterior  to  interior  zone  or  interior 
to  exterior  zone  are  distributed  from  an  origin  zone  to  a  gate  and  then 
from  the  gate  to  a  destination  zone.  This  function  requires  that  routes 
between  0-D  pairs  be  generated.  The  distribution  function  uses  travel 
time  between  0-D  pairs  as  a  parameter  in  the  equation  that  predicts 
person  trips. 

The  modal  split  function  (MODAL) ,  the  next  step  in  the  transporta¬ 
tion  planning  method,  determines  the  mode  of  travel  of  each  person,  which 
in  the  BATS  model,  is  by  civilian  vehicle,  military  vehicle,  bus,  or 
bicycle.  A  load  factor  is  determined  for  each  type  of  vehicle  operating 
from  each  zone.  This  factor  can  be  specified  as  a  function  of  O-I)  trip 
travel  time  for  buses.  Person  trips  are  then  converted  to  vehicle  trips 
using  the  load  factor  for  each  type  of  vehicle,  and  the  percentage  of 
civilian  vehicles,  military  vehicles,  buses,  or  bicycles  associated  with 
each  zone. 


The  trips  are  stored  as  the  number  of  vehicles  going  from  origin  to 
gate,  from  gate  to  destination,  and  from  internal  origin  to  internal 
destination.  At  this  point  in  the  methodological  sequence,  the  user  may 
take  advantage  of  the  fact  that  most  vehicles  on  an  air  base  travel 
through  base  gates  during  peak  hours.  The  smooth  function  (SMOOTH)  uses 
the  base  gate  counts  to  calibrate  predicted  vehicle  gate  counts  to  the 
measured  values.  Thus,  through  calibration  the  model  can  adjust  the 
predicted  number  of  vehicles  arriving  and  departing  the  base  during  any 
time  period  to  match  gate  counts.  Furthermore,  the  smf)othing  function 
shifts  trips  between  origin-gate  and  gate-destination  pairs  so  tliat  the 
predicted  number  of  vehicles  using  each  individual  gate  closely  approxi¬ 
mates  the  number  actually  counted  as  using  the  gate. 

The  final  stei)  (ASSIGN)  in  the  metliodology  assigns  vehicles  to  the 
street  network  of  the  base.  The  minimum  time  path  between  each  origin- 
gate,  gate-destination,  and  internal  origin-internal  destination  pair  is 
determined,  and  vehicle  trips  between  each  pair  are  assigned  to  the  links 
(street  segments)  in  the  route.  A  second  minimum  time  path  between  each 
zone  pair  may  be  determined  and  then  the  vehicle  trips  are  split  among 
the  two  alternate  routes  in  indirect  proportion  to  the  route  travel  times 


7 


START 


GENERATION 

FUNCTION 


DISTRIBUTION 

FUNCTION 


MODE  SPLIT 
FUNCTION 


CALIBRATION 

FUNCTION 


ASSIGNMENT 

FUNCTION 


TRIPGEN  -  PERSON  TRIPS  BY  TRIP  PURPOSE  PRODUCED  OR 
ATTRACTED  BY  EACH  ZONE  ARE  GENERATED.  FOR  HOME¬ 
WORK,  TRIP  PURPOSE,  SHIFT  COUNTS  AND  GATE  COUNTS 
INFLUENCE  THE  NUMBER  OF  TRIPS  TO  EACH  ZONE  FOR 
15-MINUTE  PERIODS. 


GRAVO  -  TRIPS  BETWEEN  ORIGIN  AND  DESTINATION 
ZONES  ARE  DISTRIBUTED.  AN  INTERMEDIATE  DIS¬ 
TRIBUTION  IS  MADE  BETWEEN  ORIGIN  ZONES  TO 
GATES,  AND  GATES  TO  DESTINATION  ZONES. 


MODAL  -  PERSON  TRIPS  ARE  CONVERTED  TO  VEHICLE 
TRIPS  USING  VEHICLE  AND  BUS  LOAD  FACTORS,  AND 
THE  NUMBER  OF  VEHICLES  OF  EACH  TYPE 


SMOOTH  -  CALIBRATION  COEFFICIENTS  ARE  DETERMINED  ^ 
BASED  ON  THE  GATE  COUNTS  DURING  THE  SIMULATED 
TIME  PERIOD.  THE  COEFFICIENTS  ARE  APPLIED  TO  ALL 
VEHICLE  TRIPS  IF  REQUESTED,  AND  MATRIX  MULTIPLICA¬ 
TION  IS  PERFORMED  TO  ROUTE  THE  APPROPRIATE  NUMBER 
OF  VEHICLES  THROUGH  APPROPRIATE  GATES.  < 


ASSIGN  -  VEHICLES  ARE  ASSIGNED  TO  LINKS  OF  THE 
NETWORK,  USING  2  MINIMUM  TIME  PATHS  BETWEEN 
EACH  ZONE-ZONE  OR  ZONE-GATE  PAIR, 


FIGURE  2.  TRANSPORTATION  PLANNING  FLOW  DIAGRAM 
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'Iln'  rouL('  travel  limes  /ire  initi/illy  l)/ist'(i  on  liit-  s|)eccl  limit  on 
e/ieli  link  in  Llie  network.  After  tiie  firj/t  iteration  of  tiie  fiA'I'S  model 
tlie  travel  times  depend  on  the  traffie  flow,  which  was  predicted  on  a 
previous  iteration.  In  this  way,  the  model  is  sensitive  to  congestion 
occurring  in  the  network,  and  routes  vehicles  away  from  areas  of  conges¬ 
tion  found  during  a  previous  Iteration. 

liATS  provides  inimts  throiigli  ASSKiN  to  ;in  <i  i  r  (piaiity  .assessment 
model.  Ttiese  inputs  t.ake  the  form  of  the  number  of  vehicles  of  six  types 
on  each  link  and  the  number  of  hot  tr.ansient,  hot  st.art,  or  cold  sl.art 
vehicles  on  e.ach  link. 

Results  of  the  five  sequential  steps  of  the  t  r.ansportat  ion  pi. inning 
methodology  are  counts  of  the  number  of  vehicles  on  each  link,  a  count  of 
the  turning  movements  made  from  each  link,  a  count  of  the  type  of  vehicle 
for  both  civilian  and  military  vehicles,  and  the  number  of  hot  transient, 
hot  start,  or  cold  start  vehicles  on  each  link.  For  the  summations  in 
the  equations  below,  the  limits  are  those  that  were  used  in  the  BATS 
computer  program. 


2.4.2  The  Trip  Generation  Function 

The  trip  generation  function  allows  the  user  to  generate  trips  from 
or  to  a  zone  in  four  different  ways  using  four  sequential  steps. 

Figure  3  shows  the  necessary  steps  and  the  tests  made  before  performing 
each  step.  The  first  step  generates  NTO,  the  trip  productions,  from  a 
zone  Z,  using  VAR,  the  demographic  variables  associated  with  zone  Z  and 
COEFO,  the  trip  productions  coefficients  associated  with  trip  purjjose  K 
and  denuigraphlc  variable  .1. 


10 

NT0(Z,K)  =  VAR(J,Z)  x  COEFO (J,K) 
J=1 


The  trip  attractions,  NTD,  to  each  zone  Z  are  similarly  generated  using 
VAR  and  COEFD,  the  trip  attraction  coefficients  associated  with  trip 
purpose  K  and  demographic  variable  J. 

10 

NTD(Z,K)  =  23  VAR(J,Z)  x  C0EFD(J,K) 

J=1 


The  second  step  in  the  trip  generation  function  is  to  determine 
default  productions  or  attractions,  if  they  are  called  for.  Default 
values  are  determined  by  using  the  land  use  designation  of  e;ich  zone 
<and  the  PLUALU  array  to  generate  trips  between  the  v.arious  land  use 
tyjies.  Each  clement  of  the  Pl.UAl.U  arr.'iy  is  thi'  numlwr  of  tot/il  trips, 
for  other  than  Home-to-Work  or  Military-Vehicle  trij)  purposes,  which  /ire 
generated  between  land  use  types.  Ttie  rows  of  the  arr.iy  .ire  Productions 
by  Land  Use  and  the  columns  are  Attractions  by  band  Use.  The  PLUAbU 
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DETERMINE  TRIP  PRODUCTIONS,  NTO,  AND  ATTRACTIONS.  NTD, 
FOR  EACH  ZONE  Z  AND  EACH  TRIP  PURPOSE  K 
all  J 

NTOIZ.K)  •  1  COEFO(J.K)  •  VARU.ZI 
all  J 

NTD(Z,KI  -  Z  COEFDIJ.K)  •  VARIJ.Z) 

NTRIPC  7 

ALUA(LI.LU)  -iZ  NTDIZ.K)  •  PLUALUILI.LUI/  Z  PLUA2U(L2.LUI 
K*2  L2-1 


No 


LU  -  LANDOd.K) 

NTOIZ.K)  -  ALUAIIZLU.LU)  •  fIVARU.Z)) 

7 

NTDIZ.K)  -  Z  ALUA(L2,LU)  •  f(VAR(J,ZII 
L2*1 


No 


COMPUTE  FLEX  TIME  POOL  SFTO  AND  SFTD 
NZONES  4 


SFTO  -  Z 

NTOIZ.II  - 

z 

SHFPCTIZ.I.I) 

Z-NEXT1 

1-1 

NZONES 

4 

SFTD  -  Z 

NTDIZ.I)  - 

Z 

SHFPCT(Z,2,I) 

Z=NEXT1 

1-1 

FIGURE  3.  FLOW  CHART  FOR  TRIP  GENERATION  FUNCTION 


JiUHU.HJI..  '  JWimnii  ijpj 


.irray  was  generated  from  an  origin-destination  study,  and  time-of-day 
tr.iv«'l  patterns  stiidic-d  at  Williams  AFI5.  An  Array  of  I, and  Use  Attractions 
lor  a  particular  time  period,  AI.IIA,  is  made  up  I  rom  the  I'l.llAl.ll  array  and 
the  N’l'lt  and  N'il)  arr.iys.  'I'lie  coliiittii  elemeni;;,  l.l,  ol  l  he  AI.IIA  .irr.iy  .ife 
filled  in  using  the  following  equation 

NTRIPC  7 

ALUA(L1,LU)  =  23  NTD(Z,K)  *  PLUALU(L1,LU)  /  PLUALU (L2 ,LU) 

K=2  /  L2=1 


where  LU  is  the  land  use  associated  with  zone  Z,  LI  goes  from  1  to  7 , 
and  NTRIPC  is  the  number  of  civilian  vehicle  trip  purposes.  The  ALUA 
array  stores  the  land  use  trips  specified  by  the  input  data  through 
the  VAR(J,Z),  and  COEF(J,K)  variables,  and  those  land  use  trips  inferred 
from  the  values  in  the  PLUALU  array. 

In  this  array  the  attractions  to  each  land  use  are  given  precedence 
over  the  productions.  Some  columns  of  the  ALUA  array  are  not  initial¬ 
ized  if  trips  have  not  been  specified  for  a  particular  land  use  on  a 
Card  Type  6.  When  this  occurs,  the  NTO  array  is  used  to  complete  all 
rows  with  land  use,  LI,  associated  with  each  Z.  The  remaining  uninitial¬ 
ized  elements  are  then  filled  in  proportion  to  the  corresponding  PLUAi.lJ 
array  elements.  The  ALUA  in  its  final  form  is  an  array  oi’  trip  produc¬ 
tions  and  attractions  associated  with  each  land  use  type. 

Default  generation  of  trips  takes  place  when  a  zero  COEFO  or  COEFD 
is  input  by  the  user  for  some  trip  purpose.  The  NTD  array  can  be  filled 
for  trip  purpose  K  for  each  zone  Z  using  its  land  use  LU  to  define  the 
total  number  of  trips  made  for  this  trip  purpose. 
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NTD(Z,K)  =  ALUA(L2,LU)  x  f(VAR(J,Z)) 
L2=l 


As  an  example,  1 f  K  is  the  shopping  trip  purpose,  then  all  zones  Z 
whose  land  use  is  shopping  (LU=3)  will  have  a  value  stored  in  NTD(Z,K). 
This  value  will  be  a  proportion  of  the  total  trips  attracted  to  all 
zones  with  land  use  =  shopping.  Because  no  particular  demographic 
variable  can  be  used  to  generate  these  attractions,  a  function  of  all 
demographic  variables  associated  with  all  zones  with  l,U=l  is  used  to 
generate  the  fraction  of  shopping  trips  going  to  a  particular  zone  Z. 

If  COEFD  is  not  zero  then; 

/all  Z 

f(VAR(J,Z))  =  COEFD(J,K)  x  VAR(J,Z)/  VAR(J,Z) 
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The  NTO  array  can  be  filled  for  trip  purpose  K,  for  each  zone,  Z,  using 
the  zones'  land  use,  IZLU,  to  determine  which  zones  produce  trips.  In 
this  case  land  use  LU  and  trip  purpose  K  are  synonymous;  that  is,  the  trip 
purpose  must  be  for  shopping,  service  or  other  land  use,  LU. 

NTO(Z,K)  =  ALUA(IZLU,LU)  x  f(VAR(J,Z)) 

1 hf  sum  of  NTO  will  therefore  equal  the  sum  of  NTD  for  any  trip  purpose  K, 
where  K=2,  l.  .  .NTRIPC. 

file  next  step  In  trip  generation  is  the  inclusion  of  shift  counts 
in  th<-  N'I'O,  Nil)  arrays  if  tiiese  are  specified  by  the  input  data.  So 
f;ir,  tiie  trip  generation  function  has  been  dealing  with  hourly  (or  perhaps 
8-hour  or  daily)  trips.  However,  to  model  congestion  during  the  peak 
periods  it  is  necessary  to  consider  15-minute  vehicle  loadings  of  the 
network  (of  course,  10- ,  5-  or  1-minute  loadings  of  the  network  would 
provide  capability  for  even  more  accurate  predictions  of  congestion,  but 
such  refinements  would  Increase  the  complications  of  the  modeling  process 
enormously) . 

A  simple  means  of  predicting  15-minute  loadings  would  be  to  make  a 
15-mlnute  iteration  of  the  BATS  model  (as  done  for  hourly  time  periods). 

The  problem  with  this  method  is  best  Illustrated  by  an  example  as  follows: 
if  zone  1  employs  100  workers  and  zone  2  employs  100  workers,  then  a 
15-minute  iteration  would  predict  an  equal  number  of  workers  going  to  both 
zones  during  any  15-minute  time  period  because  both  have  equal  demographic 
variables.  However,  75  persons  might  go  to  zone  1  and  25  to  zone  2  during 
a  15-minute  time  period.  In  order  to  represent  this  situation,  15-mlnute 
shift  counts  may  be  input  to  the  computer  and  these  will  override  the 
demographic  variable  x  coefficient  generation  of  trips  to  each  zone.  Thus, 
it  is  possible  to  input  75  employees  for  zone  1  and  25  employees  for 
zone  2.  The  shift  counts  are  input  and  stored  in  the  array  SHFPCT(Z,M, 1) 
for  each  zone  Z,  for  attractions  or  productions  M,  and  for  each  15-mlnute 
period  of  an  hour  1. 

The  final  step  is  to  use  gate  counts  to  predict  peak  15-minute 
traffic  for  other  than  shift  trips,  which  are  defined  as  FLEX  TIME 
(flexible  arrival/departure  time)  trips.  Persons  who  make  flex  time 
trips  choose  when  they  will  arrive  at  work.  A  "flex  time  pool"  is 
defined  to  include  all  non-shift  time  trips,  and  use  15-minute  gate 
counts  to  allocate  the  flex  time  pool  to  those  zones  associated  with 
gates  that  are  most  heavily  used  during  each  15  minutes. 

We  compute  GAFRO  and  GAFRD,  the  fraction  of  hourly  trips  to  a  zone 
Z  produced  or  attracted  during  a  15-mlnute  time  period: 

CAFRO  .  ^Qcm:i5iqj,x)_ 

Aft 

GC0UHT(1,G) 


GAFRD  = 


^  GCNT15(G.2,I) 
AG 

53  GC0UNT(2,G) 


where : 


p^Q  —  aXi  gates  C  associated  with  zone  Z 
GGNT15(G,J,1)  is  the  15-mlnute  gate  count  exiting 
(J=l)  through  gate  G,  during  time  period  1 
GCOUNT(J,G)  is  the  hourly  gate  count  exiting 
(J=l)  through  gate  G,  during  time  period  I. 

mtc  ttat  when  there  U  only  one  gate,  CAFRO  and  GAFRD  apply  to  all  r»ea. 
When  there  are  two  gates,  GAFRO  could  have  three  values  one  tor  zwes 
associated  with  gate  1,  one  tor  tones  associated  with  gate  2  and  one 
zones  associated  with  gates  1  and  2. 

We  next  compute  FRACO  and  FRACD,  the  fraction  of  the  hourly  flex-time 
pool  productions  or  attractions  coming  from  or  going  to  a  zone  Z. 


FRACO  = 


FRACD  = 


(nTO(Z.I)  -  E  SHFPCT(Z,l,I)j 

all  IZ'/  \ 

53  /nT0(1Z,1)  -  5]  SHFPCT(IZ,1,1) j 

|nTD(Z,1)  -  E  SHFPCT(Z,2,I) J 

all  IZ  /  Jt,  A 

53  /ntd(iz,i)  -  E  SHFPCT(1Z,2,I)| 


FRACO  X  GAFRO  x  SFTO  is  the  fraction^ of  the  hourly  flex-time  pool, 
SFTO,  coming  from  zone  Z  during  a  15-minute  period. 

all  IZ  /  Jl,  \ 

where  SFTO  =  E  (nTO(IZ.I)  -  E  SHFPCT(IZ,1 ,1) j 

SFTD  =  E  |nTD(TZ,1)  -  E  SHFPCT(1Z,2,I)| 

The  number  of  trips  coming  from  any  zone  Z  during  a  15-mlnute  time-  period 
I  is  then: 
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NT0(Z,1)  =  FRACO  x  SFTO  x  GAFRO  +  SHFPCT(Z ,1 , I) 


The  number  of  trips  going  to  any  zone  Z  during  a  15-minute  time 
period  I  is: 

NTD(Z,1)  =  FRACD  x  SFTD  x  GAFRD  +  SHFPCT(Z ,2 , I) 


2.4.3  Trip  Distribution  Function 

Thus  far,  the  NTO(Z,K)  and  NTD(Z,K)  arrays  have  been  developed. 

These  arrays  store  the  person-trips  originating  from  or  destined  to  a 
zone  Z  for  the  trip  purpose  K.  The  distribution  function  determines 
the  proportion  of  the  NTO(Z,K)  trip  origins  going  to  every  other  zone 
in  the  network,  and  it  determines  the  proportion  of  the  NTD(Z,K)  trip 
destinations  coming  from  every  other  zone  in  the  network. 

One  of  the  widely  used  methods  of  distributing  trips  is  through  the 
use  of  a  "gravity"  function,  so  named  from  Newton's  second  law  of  motion, 
which  states  that  the  attraction  between  two  bodies  is  directly  propor¬ 
tional  to  the  product  of  their  masses  and  indirectly  proportional  to  the 
sf|uar€!  of  the  distance  between  them.  The  trip  distribution  function 
assumes  that  the  trip  attraction  between  two  zones  is  directly  propor¬ 
tional  to  the  product  of  the  number  of  trip  origins  and  trip  destinations 
and  indirectly  proportional  to  a  function  of  the  travel  time  between  them. 

On  an  air  base,  another  variable  must  be  considered  in  the  function — 
the  number  of  vehicles  traveling  through  each  gate  to  the  base.  More 
than  one  trip  travel  time  between  zones  is  possible;  each  trip  may  go 
through  one  or  more  gates  on  the  base.  The  problem  faced  in  defining  a 
distribution  function  for  an  air  base  was  to  use  the  gate  counts  to  get 
a  more  accurate  distribution  of  trips  and  to  route  the  correct  number  of 
trips  between  each  pair  of  zones. 

Previous  studies  of  traffic  at  the  air  bases  (Nellis,  1971; 

Tinker,  1967)  provided  data  that  could  be  used  to  determine  the  function 
of  gate  counts  and  travel  time  on  vehicle  distribution.  A  study  of  trip 
travel  patterns  at  Nellis  AFB  showed  that  trips  made  from  gate  to  employ¬ 
ment  zones  were  Indirectly  proportional  to  travel  time  cubed. 

The  Nellis  AFB  study  predicted  the  number  of  vehicles  traveling 
from  each  of  four  gates  to  each  of  twelve  work  zones.  The  number  of 
vidiicles  gf)ing  from  each  g<ate  to  each  destination  zone  was  projected 
I  rom  interviews  with  ap|)roximate)  y  one-third  of  the  vehicles  using  tiie 
gates  during  the  morning  rush  hour  and  documented  In  the  1971  Traffic 
Knglneerlng  Study. 
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Thf  Htiulli'H  at  rinkcT  AFIi,  Nfllis  AKB,  and  other  b  as<-s  showed  that 
employees  fended  to  use  a  gate  near  to  llu-lr  zone  of  employmet)!  and 
employees  tended  ncjt  to  double  back  in  their  route,  i  .  i- .  ,  drive  beyond 
an  employment  zone  while  off  base  and  back  to  the  zone  while  on  base. 

One  modification  to  the  distribution  procedure  was  suggested  by  these 
studies.  This  would  lower  attraction  to  a  zone  from  a  gate  when  trips 
within  the  base  retrace  their  trip  while  approaching  the  base;  in  other 
words,  when  vehicles  must  drive  past  their  destination  zone  in  order  to 
use  a  certain  entrance  gate. 

A  trip  distribution  function  that  is  inversely  proportional  to  off- 
base  travel  time  plus  on-base  travel  time  cubed,  and  directly  proportional 
to  gate  counts  appears  to  satisfy  the  constraints  of  the  problem.  The 
equations  used  in  the  GRAVO  subroutine  (which  performs  the  trip  distribu-- 
tion  function)  are  described  in  the  following  paragraphs.  A  flowchart 
of  the  function  is  given  in  Figure  A. 

The  two  ends  of  the  trips  are  in  exterior  zones  Z  and  interior 
zones  Zl;  that  is,  the  equations  handle  trips  that  originate  in  Z  and  go 
to  Zl  and  that  originate  in  Zl  and  go  to  Z. 

The  number  of  trips  between  Z  and  Zl  is  considered  to  be  directly 
proportional  to  the  fraction  of  the  total  destined  for  Zl.  (This 
assumes  that  employees  do  not  choose  their  housing  location  based  on 
their  zone  of  employment.) 

all  K  /all  Zl 

ZGT  =  52  NT()(Z,K)  ^  NTI)(ZI,K)/  'NTI)(Z1,K) 

all  K  /all  Zl 

ZAT  =  52  NTD(Z,K)  x  NT0(Z1,K)/  5Z  NT()(ZI,K) 


Then  ZGT  is  the  total  number  of  trips  leaving  from  zone  Z  for  zone  Zl 
and  ZAT  is  the  total  number  of  trips  going  to  zone  Z  from  zone  Zl. 

Let  TT(Z,G)  be  the  travel  time  from  zone  Z  to  gate  G.  Then  the 
fraction  of  ZGT  or  ZAT  that  uses  any  gate  is  proportional  to  FSODl  or 
FS0D2 : 


FSODl(G)  = 


GCNT15 (G,1,ITM) 
TT(Z,G)  +  TT(Z1,G)^ 


FS0D2(G)  = 


GCNT15(G,2,ITM) 
TT(Z,G)  +  TT(Z1,G)^ 


where  GCNT15(G,1,ITM)  is  the  15-mlnute  gate  count  exiting  gate  G  during, 
the  ITM  time  period,  and  GGNTl 5 (G ,2 , ITM)  is  the  15-minute  gate  count 
entering  gate  G.  So  ZGT  x  FSODl(G)  is  the  gravity  function  applied  tt> 
any  gate  for  a  pair  of  zones  Z  and  Zl. 
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INO  SET  POSITIVE 


INO  SET  NEGATIVE 


DETERMINE  FSOD(G),  THE  FRACTION  OF  TRIPS  GOING  BETWEEN 
EXTERNAL  ZONE  Z  TO  INTERNAL  ZONE  Z1  THROUGH  GATE  G 

GCNT15(G,1,ITM) 

TTIZ.G)  +  TT(Z1,G)3 

GCNT15(G.2.ITMI 

TTIZ.G)  +  TTIZI.G)^ 

~  ■  '  ■  ■  .1  ■  ,  , 

DETERMINE  OV(Z,G)  AND  VDIZ.GI,  THE  NUMBER  OF  TRIPS 
GOING  FROM  ZONE  Z  TO  GATE  G  OR  FROM  GATE  G  TO  ZONE  Z 

OVIZ.GI  =  FSOD2(G)  •  ZGT/  1'  FSOD2(GI 
OVIZI.GI  =  FSODKG)  *  ZAT/  X  FSODKG) 

VDIZ.GI  =  FSODKG)  *  ZAT/  X  FSODKG) 

VDIZI.G)  =  FSOD2(G)  *  ZGT/  X  FSOD2(G) 


tor  all  trip  purposes  K 

I 

for  all  internal  zones  Z1 


DETERMINE  INTERNAL  TO  INTERNAL  ZONE  TRIP 
PRODUCTIONS  (NG1  IS  THE  NUMBER  OF  GATES  +  1) 

(all  Z1  alt  Z 

_X  I^OIZI.K)  -  X  NTD(Z.K) 

all  Z1 

X  NTOIZt.K) 


H 


FIGURE  4.  FLOW  CHART  FOR  TRIP  DISTRIBUTION  FUNCTION 
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FIGURE  4. 


FL^W  CHART  FOR  TRIP  DISTRIBUTION  FUNCTION  (Concludud) 
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The  function  still  must  be  normalized  so  that  tlie  sum  of  the 
ZGT  X  FSODl(G)  will  equal  Z(;T.  In  ])errorming  the  normalization  tiie 
model  computes  the  number  going  from  zone  Z  to  each  gate  G,  OV(Z,G),  the 
number  going  from  zone  Z1  to  each  gate  G,  0V(Z1,G),  the  number  going  to 
zone  Z  from  gate  G,  VD(Z,G),  and  the  number  going  to  zone  Z1  from  gate  G, 
VD(Z1,G) . 


OV(Z,G)  =  ZGT  X  FS0D2(G)  ^ 

/all  G 

/  E 

FS0D2(G) 

0V(Z1,G)  =  ZAT  X  FSODl(G)  j 

'  E 

FSODl(G) 

VD(Z,G)  =  ZAT  X  FSODl(G) j 

/all  G 

'  E 

FSODl(G) 

VD(Z1,G)  =  ZGT  X  FS0D2(G)  j 

1  all  G 

1  E 

FS0D2 (G) 

The  result  of  the  distribution  function  GRAVO  is  then  the  OV  and 
VD  arrays  which  store  an  estimate  of  the  number  of  person  trips  from 
any  zone  to  any  gate  or  from  any  gate  to  any  zone.  Because  of  dealing 
with  person  trips,  the  results  are  an  estimate.  Later  in  the  MODAL  and 
SMOOTH  subroutines,  the  OV  and  VD  arrays  are  further  adjusted  so  that 
they  more  closely  approximate  the  gate  counts. 

all  Z 

GCNT15(G,1,ITM)  =  OV(Z,G) 


and 


all  Z1 

GCNT15(G,2,1TM)  =  ^  VD(Z1,G) 


The  prediction  of  trips  from  internal  zone  to  internal  zone  is  stored 
in  the  0V(Z1,NG1)  and  VD(Z1,NG1)  arrays,  where  NGl  is  the  number  of 
gates  plus  1.  The  internal  trips  are  those  originating  from  internal 
zones  that  are  not  attracted  to  external  zones  and  those  destined  to 
internal  zones  that  df)  not  originate  in  external 

all  Z1  all  K  /all  Z1 

Y,  0V(Z1,NG1)  =  2]  \  I]  NT0(Z1,K) 

all  Z1  all  K  /all  Z1 

Y,  VD(Z1,NG1)  =  £  y  2]  NTD(Z1,K) 

The  internal  trips  will  be  assigned  later  by  trip  purpose  by  the 
assignment  function.  Thus,  the  NTO  and  NTD  arrays  are  updated  to  repre¬ 
sent  internal  trips  by  trip  p\irpose.  The  equations  that  redefine  NTO 
and  NTD  are: 


zones . 


all  Z 


Z  \ 

23  NTD(Z,K)| 

all  Z  \ 

-  X)  NTO(Z,K)j 
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all  Z1  all  Z 

_  2^  NT0(Z1,K)  -  ^  NTD(Z.K)  x  NTO(Zl.K) 

^  ^  all  Z1 

NT0(Z1,K) 

all  Z1  all  Z 

NTD (Z1 , K)  =  .  .X..  _iLTIjL(.ZJ^j<)_ji_  .  N:[0 

NTO(Zl,K) 

2.4.4  Mode  Split  Function 

The  mode  split  function  is  the  third  major  step  in  the  transportation 
planning  methodology,  and  converts  person  trips  to  vehicles.  Vehicle  load 
factors  are  used  to  determine  the  number  of  autos  and  trucks  used  to  make 
trips.  Bus  and  bicycle  traffic  is  handled  relatively  easily;  bicycles 
are  considered  to  have  a  load  factor  of  1,  and  buses  have  a  load  factor 
that  is  specified  by  input  data.  Figure  5  is  a  flowchart  of  the  mode 
split  function. 

The  number  of  bus  passengers  is  subtracted  from  the  OV,VD  demands 
as  the  first  step  in  this  subroutine.  The  bus  load  factor  is  given  by 
VLFM(7)  and  the  number  of  buses  associated  with  any  zone  Z  is  V1TPM(7,Z). 
Thus,  the  number  of  persons  using  the  bus  when  leaving  any  zone  is: 


FI  =  VLFM(7)  X  VTYPM(7,Z) 


It  is  possible  to  let  the  computer  determine  the  bus  load  factor  based 
on  the  number  of  buses  servicing  the  zone  and  the  number  of  trips  from 
the  zone.  If  VLFfI(7)  is  zero,  then  the  number  of  bus  passengers  in  an 
hour  is: 


FI  =  VTYPM(7,Z)  X  (0V(Z,NG1)  +  VD(Z ,NG1) ) /36 


The  constant  36-person  trips  bus  trip  was  derived  to  fit  the  general 
Information  available  on  bus  loadings  on  air  bases. 

The  load  factor  for  each  type  of  vehicle,  I,  is  given  by  VI.F(T)  for 
civilian  vetiicle  trips  or  VLFM(I)  for  military  vehicle  trips.  VFUi'YP  ( I  ,  Z) 
givc.'s  the  percent  of  vehicles  of  type  I  associated  with  zone  Z  and 
VHTYPM(I,Z)  gives  the  percent  of  military  vehicles  of  typi'  I; 
all  I 

2  VEHTYP(I,Z)  can  be  less  than  100  and  this  difference  represents 
the  percentage  of  nonmotorized  trips  made  from  or  to  the  zone. 

To  compute  vehicle  trips  for  zone  Z,  we  first  determine  SPPHl  and 
SPPH2,  the  number  of  persons  occupying  100  vehicles  (bicycle  riders  are 
included) . 
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COMPUTE  PERSONS  USING  BUSES  IN  EACH  ZONE 
FI  -  VLFM(7)  •  VTYPM(7.ZI 


T 


FIGURE  5.  FLOW  CHART  FOR  MODE  SPLIT  FUNCTION 


all  I 

all  I 

SPPHl  =  53 

VLF(l)  X 

VEHTYP(I,Z)  +  100  - 

53  VEHTYPd, Z) 

all  I 

all  1 

SPPH2  = 

VLFM(I)  X 

VHTYPM(I,Z)  +  100 

-  53  VHTYPMd,Z) 

Then  person  trips  are  converted  to  vehicle  trips  in  the  OV,  VD, 
NTO  and  NTD  arrays: 


all  I 

NTO(Z,K)  =  NTO(Z,K)  x  VEHTYP(l,Z)/SPPHi 

all  I 

NTD(Z,K)  =  NTD(Z,K)  x  VEHTYP(I, Z) /SPPHl 

,  all  I 

OV(Z,G)  =  OV(Z,G)  X  |FRAM1L(Z)  x  VHTYPM(I ,Z) /SPPH2 

all  I  V 

+  (1  -  FRAMIL(Z))  X  53  VEHTYP(I,Z)/SPPHll 

y  all  I 

VD(Z,G)  =  VD(Z,G)  X  fFRAMIL(Z)  x  V11TYPM(I , Z) /SPPH2 

V  all  I 

+  (1-FRAMIL(Z))|  X  53  VEHTYPd  ,Z) /SPPHl 

where  FRAMIL(Z)  is  the  fraction  of  military  vehicles  making  trips  through 
the  gates  from  or  to  zone  Z. 


2.4.5  Calibration  Using  Gate  Counts 

The  unique  fact  about  an  air  base  that  can  contribute  to  the  accuracy 
ol  the  model  is  that  external  traffic  must  travel  through  gates  to  get 
on-base.  The  gate  counts  for  each  base  gate  act  the  same  as  "cordon 
line"  counts,  which  are  typically  used  in  a  transportation  planning 
process  to  calibrate  a  model.  In  calibration,  modeled  results  are  cor¬ 
rected  to  more  closely  approximate  results  obtained  from  field  measure¬ 
ment.  In  BATS  this  function  serves  to  account  for  those  trips  that  were 
not  modeled  using  the  data  inputs.  A  flowchart  of  the  principal  sub¬ 
routine  SMOOTH  is  shown  in  Figure  6. 

From  previous  routines  or  functions,  the  OV(Z,G),  VD(Z,G)  arrays 
that  store  the  vehicle  trips  going  from  zone  Z  to  gate  G  and  from  gate 
G  to  zone  Z,  respectively,  have  been  established.  Defining  all  Z  as 
exterior  zones  and  Z1  as  Interior  zones,  calibration  coefficients  are 
determined  as  follows: 


all  Z  all  G  /all  G 

FEXGEN(ITM)  =  2}  53  OV(Z,G)/  «CNT15  (G  ,2 ,  ITM) 

where  FEXGEN  is  the  fraction  of  external  generations  created  by  the  model 
and  traveling  through  the  gates  and  GCNT15(G ,2 ,1TM)  are  the  15-minute  (or 
hourly)  counts  of  the  vehicles  entering  gate  G  during  time  period  ITO 
(1TM*1  for  hourly  iterations). 
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^  SMOOTH  ^ 


^  For  all  gates  G 


FIGURE  6.  FLOW  CHART  FOR  CALIBRATION  FUNCTION 
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FACTOR  THE  OV  AND  VD  ARRAYS  SO  THE  COLUMNS 
SUM  TO  THE  GATE  COUNTS. 

all  Z 

OV(Z,G)  »  OV(Z,Gl  *  GCNT16(G,2,ITM)/  I  OV(Z,GI 

all  Z1 

OV(Z1,G)  »  OV(Z1,G)  •  GCNT15(G,1.ITMI  /  I  OV(Z1,GI 

all  Z 

VDIZ.GI  '  VD(Z,G)  *  GCNT15(G.1.ITMI  /  Z  VD(Z.G) 

all  Z 

VO(Z1,G)  =  VD(Z1,G)  •  GCNT15(G,2,ITMI  !  2  VD(ZI.G) 


FIGURE  6.  FLOW  CHART  FOR  CALIBRATION  FUNCTION  (Continued) 
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similarly,  FEXATT,  the  fraction  of  the  total  external  attractions 
created  by  the  model,  FINGEN,  the  fraction  of  the  total  Internal  genera¬ 
tions  created  by  the  model,  and  FINATT,  the  fraction  of  the  total  internal 
attractions  created  by  the  model  are  computed. 

all  Z  all  G  /all  G 

FEXATT(ITM)  =22  VD(Z,G)/  2  GCNT15 (G , 1 , ITM) 
all  Z1  all  G  /all  G 

FINGEN(ITM)  =2  2  0V(Z1,G)/  ^  GCNT15 (G , 1 , ITM) 

all  Z1  all  G  /all  G 

FINATT(ITM)  =22  VD(Z1,G)/  2  GCNT15 (G,2 ,ITM) 

If  the  user  chooses  to  calibrate  the  model  to  the  gate  counts,  then 
the  elements  of  the  OV  and  VD  arrays  are  updated  by  applying  the  FEXGEN, 
FEXATT,  FINGEN,  and  FINATT  factors  as  follows: 

OV(Z,G)  =  OV(Z,G)/FEXGEN(ITM) 

VD(Z,G)  =  VD(Z,G)/FEXATT(ITM) 

0V(Z1,G)  =  0V(Z1,G)/FINGEN(ITM) 

VD(Z1,G)  =  VD(Z1,G)/ FINATT (ITM) 

The  user  selects  the  option  making  external-to-lnternal  and  intemal- 
to-external  trip  productions  and  attractions  total  to  the  gate  counts. 

The  user  also  has  the  option  of  making  external  productions  total  to 
internal  attractions,  or  of  making  internal  productions  total  to  external 
attractions.  The  user  may  also  input  the  FEXGEN,  FEXATT,  FINGEN,  and 
FINATT  values  when  gate  counts  are  not  available  for  a  predictive  run. 

The  GRAVO  subroutine  predicts  the  number  of  person  trips  through  each 
gate,  using  the  gate  counts  at  each  gate  as  a  weighting  factor.  The  total 
number  of  entering  or  exiting  vehicles  now  sums  to  the  gate  counts  if  the 
calibration  factors  have  been  applied.  Further  calibration  adjusts  the 
OV  and  VD  arrays  to  gate  counts  so  that  the  sum  of  the  OV  going  through 
any  gate  G  more  nearly  equals  the  count  at  gate  G.  This  step  is  accom¬ 
plished  using  a  matrix  multiplication  scheme.  Travel  time  is  not  con¬ 
sidered  to  be  a  factor  any  longer,  and  the  number  using  each  gate  is  most 
important . 

The  matrix  multiplication  is  a  two-step  procedure:  (1)  factor  all 
elements  so  that  the  column  sums  for  external  zone  productions  add  up  to 
the  entering  gate  counts,  the  column  sums  for  Internal  zone  productions 
add  up  to  the  exiting  gate  counts,  the  column  sums  for  external  attrac¬ 
tions  add  up  to  exiting  gate  counts,  and  the  column  sums  for  internal 
attractions  add  up  to  entering  gate  counts;  (2)  factor  all  elements  so 
that  row  sums  for  each  zone  equal  the  initial  row  sums  for  each  zone. 

Thus,  the  initial  row  sums  of  productions  from  each  zone  are  saved: 

all  G 

ZGENR(Z)  =  2  OV(Z,G) 


The  elements  are  factored  so  that  the  columns  sum  to  the  gate  counts: 


/all  Z 


OV(Z,G)  =  OV(Z,G)  X  GCNT15(G,2,ITM)/  OV(Z,G) 


/all  Z1 


0V(Z1,G)  =  0V(Z1,G)  X  GCNT15(G,1,ITM)/  0V(Z1,G) 


where  Z  are  external  zones  and  Z1  are  internal  zones.  The  elements  are 
factored  so  that  row  sums  for  each  zone  equal  initial  row  sums. 


OV(Z,G)  =  OV(Z,G)  X  ZGENR(Z) 


'all  G 

E 


OV(Z,G) 


Similarly  for  zone  attractions,  the  initial  row  sums  to  each  zone  are 
saved : 


all  G 

ZATTR(Z)  =  53  VD(Z,G) 

The  elements  are  factored  so  that  the  columns  sum  to  the  gate  counts: 

/all  Z 

VD(Z,G)  =  VD(Z,G)  /  GCNT15(G,1,1TM)/  VIJ(Z,G) 

/all  Z1 

VD(Z1,G)  =  VD(Z1,G)  X  GCNT15(G,2,ITM)/  VD(Z1,G). 

/all  G 

VD(Z,G)  =  VD(Z,G)  X  ZATTR(Z)/  VD(Z,G) 

Four  iterations  of  this  procedure,  using  Nellis  AFB  data,  produced 
OV  and  VD  matrices  in  which  the  column  sums  were  nearly  equal  to  the  gate 
counts.  BATS  is  programmed  to  perform  the  matrix  multiplication  four 
times. 

The  final  step  in  the  calibration  function  is  to  reroute  vehicles 
to  alternative  lots,  away  from  parking  lots  that  are  filled  to  capacity. 
The  ZGENR(Zl)  array  is  made  to  store  all  trips  produced  from  each  internal 
zone  Zl,  and  the  ZATTR(Zl)  array  is  made  to  store  all  attractions  to  each 
internal  zone  Zl  by  adding  the  intemal-to- internal  zone  trips  to  these 
arrays : 


all  K 

ZGENR(Zl)  =  ZGRNR(Zl)  +  NT0(Z1,K) 

all  K 

ZATTR(Zl)  =  ZATTK(Zl)  +  NTD(Zl,K) 
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where  NTO  and  NTD  store  the  number  of  internal  origins  going  to  internal 
destinations  and  the  number  of  internal  destinations  coming  from  inter¬ 
nal  origins,  respectively.  The  vehicles  that  use  parking  places  are  of 
type  1,  2  and  3,  so  we  generate  ZM(Zl)  and  VZ(Zl) .  the  number  of  parked 
vehicles  originating  in  zone  Z1  and  the  number  of  vehicles  wishing  t(j 
park  in  zone  Zl,  respectively. 


3  /all  1 

ZV(Zl)  =  ZGENR(Zl)  X  5^  VEHTYPd.Zl)/  ^  VEHTYP(I,Z1) 

1=1  ' 

3  /all  I 

VZ(Zl)  =  ZATTR(Zl)  X  VEHTYP(I,Z1)/  VEHTYP(I,Z1) 

1=1 

If  the  capacity  of  the  lot  is  exceeded,  then  excess  vehicles  are  routed 
to  a  parking  zone  that  has  a  link  in  common  with  the  overcapacity  zone, 
and  the  VZ  and  ZATTR  trip  attractions  arrays  are  updated  to  reflect  this 
rerouting.  The  rerouting  is  made  to  an  adjacent  zone,  or  if  none  exists 
to  the  zone  with  the  greatest  unused  capacity. 


2. A. 6  Assignment  Function 

This  part  of  the  transportation  planning  methodology  assigns  vehicles 
to  the  links  of  the  network.  Vehicle  trips  are  assigned  to  the  minimum 
path  routes  between  each  origin  and  gate,  each  gate  and  destination,  and 
each  internal-internal  zone  pair.  Vehicle  type  data,  including  civilian 
or  military  vehicles,  are  maintained  for  each  link.  Data  on  hot,  cold, 
and  stabilized  vehicles  are  also  accumulated  for  eacli  link.  When  a  truck 
route  has  been  specified,  the  vehicles  of  the  appropriate  type  are  routed 
on  the  truck  route  instead  of  the  minimum  patli  route.  The  arrays  that 
store  the  minimum  path  and  the  minimum  path  cost  (or  travel  time)  are  the 
R(L),  C(L),  RG(G,L)  and  CG(G,L). 

The  R  array  stores  the  link  previous  to  L  in  the  route  from  the  zone 
of  origin  to  link  L.  Thus,  if  R(6)  =  4,  the  route  from  the  zone  of  origin 
to  link  6  passes  over  link  A.  When  R(L)  =  0,  there  are  no  other  links  in 
the  route,  so  the  link  L  accesses  the  origin  zone.  The  RG(G,L)  array  is 
like  the  R(L)  array,  only  it  stores  the  link  previous  to  L  in  the  route 
from  gate  G  to  link  L.  The  cost  array  C(L)  stores  the  cost  (or  travel 
time)  from  the  origin  zone  to  link  L,  and  CG(G,L)  stores  the  cost  (or 
travel  time)  from  gate  G  to  link  L.  Figure  7  is  a  flowchart  of  the 
assignment  subroutine  ASSIGN. 

The  assignment  takes  place  in  four  steps: 

(1)  For  each  zone  Z  the  R(L)  and  C(L)  arrays  are  determined. 

(2)  For  each  gate  G  the  0V(Z,G)  vehicles  are  routed  onto  the  links 
of  the  path  from  zone  Z  to  gate  G,  and  the  VD(Z,G)  vehicles 
are  routed  from  the  gate  G  onto  the  links  of  the  two  alteni.iic 
paths  to  zone  Z. 
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M 


For  each  /one  Z 


'CALL  MNPATH  TO  FIND  MINIMUM  PATH  FROM  Z> 
TO  ALL  LINKS  L  OF  THE  NETWORK 
STORE  RESULTS  IN  R(L)  AND  C(L) 


- 1 - 

For  each  gate  G  ' 


FIGURE  7.  ASSIGNMENT  OF  VEHICLES  TO  NETWORK 
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} 


For  all  internal  zones  Z1 
For  each  trip  purpose  K 


FIGURE  7.  ASSIGNMENT  OF  VEHICLES  TO  NETWORK  (Concluded) 
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(3)  For  each  trip  purpose  K  the  NTO(Z,K)  vehicles  are  routed  onto 
the  two  alternate  paths  to  all  internal  zones  Zl.  The  number 
routed  from  Z  to  Zl  (OD)  is  proportional  to  the  attractions 
to  zone  Zl  for  the  Kth  trip  purpose. 

/all  Zl 

OD  =  NTO(Z,K)  X  NTD(Z1,K)/  NTD(Z1,K) 

(4)  The  trips  OD  between  zone  Z  and  Zl  are  divided  among  the  two 
alternate  paths  as  follows: 

Trips  ending  on  link  L2  =  OD  x  C(L1)/[C(L1)  +  C(L2)] 

Trips  ending  on  link  LI  =  OD  x  C(L2)/[C(L1)  +  C(L2)] 

Steps  1  through  4  are  repeated  for  each  zone  Z  at  which  point  all  vehicles 
have  been  assigned  to  the  network. 

The  results  are  stored  In  the  COUNT(M,L)  array  where  L  is  a  link 
number  and  M  is  the  index  for  through  (M=l) ,  right  (M=2) ,  left  (M=3) , 
terminating  (M=4) ,  vehicle  type  1-6  (M=5-10) ,  military  vehicle  type  1-6 
(M=ll-16) ,  number  of  cold  starts  (M=17) ,  and  number  of  hot  starts  (M=18) . 

The  ASSIGN  subroutine  performs  the  assignment  function,  and  it  con¬ 
sumes  far  more  computer  time  than  any  other  subroutine  in  the  BATS  com¬ 
puter  model.  However,  the  user  may  minimize  the  computer  run  time  using 
the  I0PT(5)  option.  This  option  allows  the  subroutine  to  follow  only  one 
path  from  origin  to  destination  zones  and  to  route  small  numbers  of 
vehicles  traveling  from  Z  to  Zl  to  alternative  zones  so  that  fewer  paths 
are  followed.  Thus,  from  the  dimensions  of  the  arrays,  there  are: 

50  origin  zones  x  iQ  gates  =  500  possible  paths  + 

50  destination  zones  x  lo  gates  =  500  possible  paths  + 

44  internal  zones  x  43  internal  destinations  =  1892  possible  paths  + 

50  destination  zones  x  10  gates  =  500  possible  alternate  paths  + 

44  internal  zones  x  43  internal  destinations  =  1892  possible  alternate  paths. 

There  can  be  approximately  4500  paths  followed  through  the  net¬ 
work  during  the  trip  assignment  process,  and  if  following  each  path 
takes  10  msec  it  would  be  possible  to  use  45  seconds  of  computer  time 
just  assigning  vehicles  to  the  network  for  one  time  period.  A  run  of  12 
hourly  time  periods  with  one  morning  and  one  evening  rush  hour  simulated 
by  four  15-minute  time  periods  could  use  more  than  13  minutes  of  computer 
run  time. 
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2 . 5  Traffic  Flow  Analysis 
2.5.1  General  Methods 

General  traffic  flow  analysis  covers  such  subjects  as  network  theory, 
traffic  assignment,  queueing  theory,  interrupted  and  uninterrupted  flow 
theory,  surveillance,  and  control.  The  analysis  of  traffic  in  and  about 
an  indirect  source  such  as  an  air  base  principally  Involves  traffic 
assignment,  queueing  theory,  and  street  and  parking  lot  flow.  Traffic 
analysis  at  an  air  base  may  be  based  on  simplistic  methods,  which 
generally  use  hand  calculations,  or  on  sophisticated  methods,  which 
generally  Involve  computer  modeling.  Hand  calculations  used  in  traffic 
flow  analysis  suffer  from  being  cumbersome  and  insensitive  to  some  major 
parameters  of  a  typical  traffic  network. 

Manual  analytical  procedures  for  vehicle  behavior  in  and  around 
Indirect  sources,  such  as  an  air  base  with  most  pollution  generated  by 
aircraft  or  motor  vehicles  (EPA  Guidelines,  1978;  Thayer  and  Axetell, 
1973),  recognize  the  need  to  model  individual  segments  of  a  vehicle's 
journey,  but  the  modeling  of  intersections,  parking  lot  running  time, 
and  vehicle  routing  into  and  out  of  a  parking  lot  are  very  simplified. 

The  simplification  is  justifiable  to  keep  hand  computations  to  a  minimum, 
but  such  simplification  may  lead  to  significant  inaccuracies  in  modeling 
the  unusual  circumstances  that  are  typical  of  worst-case  conditions. 

For  example,  assune  a  fairly  congested  exit  gate  that  leads  to  a  signal¬ 
ized  intersection  operating  at  near  capacity  with  a  cycle  length  of  1.5 
minutes  and  a  red  phase-to-cycle-length  ratio  of  0.5.  A  vehicle-actuated 
signal  might  operate  at  a  1.5-minute  cycle  length  but  could  allocate  more 
time  to  the  green  phase,  thus  effectively  lowering  the  volume-to-capacity 
ratio  and  the  travel  time  through  the  gate.  The  hand  computation  tech¬ 
niques,  however,  contain  no  provision  for  determining  phase  and  cycle 
lengths  of  actuated  signals. 

From  previous  experience  in  modeling  traffic  networks  using  a  large- 
scale,  complex  computerized  model  (Sandys,  1971;  Sandys  et  al.,  1975), 
the  parameters  that  are  the  most  important  in  prediction  of  flow  on  a 
netwoik  of  streets  have  been  isolated.  Of  major  importance  are  inter¬ 
section  capacity,  number  of  lanes,  and  number  of  left-turning  vehicles; 
of  If.'ss  Importance  are  street  capacity  and  parking  characteristics.  On 
the  basis  of  the  major  parameters,  a  simple  street  network  flow  model 
has  been  developed  that  automatically  computes  travel  times  based  on 
trip  demand,  intersection  capacities,  optimum  phase  lengths  for 
actuated  signals,  and  parking  lot  travel  times. 

The  traffic  flow  submodel  is  rerun  for  each  time  period  that  is 
specified  (see  Figure  8).  A  time  period  is  specified  as  an  input  param¬ 
eter  and  a  situation  is  modeled  during  this  time.  Accumulations  of 
vehicles  in  parking  lots  and  queues  on  streets  are  retained  and  during 
subsequent  time  periods  these  act  as  initialization  data  for  the  subse¬ 
quent  iteration. 
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CALL  INSECU  TO 
COMPUTE 

INTERSECTION  DELAY 
AND  QUEUE  LENGTH 


<  0 


CALL  INSEC  TO  COMPUTE  INTERSECTION 
delays,  queue  LENGTHS  AND  NUMBER  OF 
STOPS  DURING  THE  ITERATION  TIME  PERIOD 


CALL  PARKING  TO  COMPUTE  THE  RUNNING  TIM^ 
V  OF  VEHICLES  USING  EACH  PARKING  ZONE  y 


- - 

SAVE  AVERAGE  LINK  SPEEDS,  DELAYS,  QUEUE 
LENGTHS,  VOLUMES  AND  NUMBER  OF  STOPS, 

AND  SAVE  PARKING  ZONE  RUNNING  TIMES  AND 
VOLUMES  FOR  USE  BY  EMISSION  AND 
DIFFUSION  MODELS 

RETURN 


FIGURE  8.  TRAFFIC  FLOW  ANALYSIS 
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2,5,2  Intersections 


Traffic  is  assigned  to  the  altentate  routes  on  the  basis  of  time 
delay  and  capacity;  the  resulting  travel  times  and  other  measures  of 
effectiveness,  such  as  stops,  delay,  queue  length,  and  flow,  arc  then 
calculated.  The  delays  at  each  intersection  are  calculated  by  a  sub¬ 
routine  that  models  an  intersection  on  the  basis  of  traffic  flow  and 
physical  characteristics.  The  intersection  is  modeled  in  a  simple 
geometrical  model  (Figure  9)  to  permit  many  intersections  and  demand 
patterns  to  be  simulated  in  a  short  time.  The  outputs  of  tlu-  subroutine 
are  average  vehicle  delay  and  average  queue  length  for  eacii  apptoacli  and 
turning  movement  at  the  intersection. 


2, 5, 2,1  Signalized  Interactions 

The  method  of  computing  delay  at  a  signalized  intersection  is  based 
on  an  approximate  method  of  computing  delays  and  queues  (Newell,  1965) 
in  which  the  discrete  nature  of  the  cars  is  disregarded  in  favor  of  con¬ 
sidering  traffic  as  a  continuous  fluid,  which  arrives  at  a  uniform  rate 
(V) ,  is  dammed  for  a  time  (R) ,  and  is  then  released  at  a  rate  (S)  until 
the  dam  is  empty.  Traffic  thereafter  moves  out  of  the  intersection  at 
the  arrival  rate  (V)  provided  the  green  time  (G)  is  long  enough. 


VEH 


l^R  - -  0  - ^ 

FIGURE  9.  SIMPLE  MODEL  OF  INTERSECTION  DELAY 
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In  Figure  9,  the  initial  form  of  the  subroutine,  VEH,  is  the  vehicle 
arrivals.  The  total  delays  for  all  vehicles  during  any  time  interval  are 
represented  by  the  interval  bounded  by  the  arrival  and  departure  flow 
rates.  The  average  delay  (W)  per  vehicle  is: 

W  =  S*R*R/(2*(R  +  G)*(S  -  V)) 

Actually,  vehicles  rarely  arrive  at  a  uniform  rate  but  tend  to 
arrive  in  random  small  groups.  In  the  model,  a  Poisson  distribution  has 
been  assigned  to  this  characteristic,  and,  when  such  arrivals  occur, 
queues  of  vehicles  form  waiting  to  be  serviced  by  the  intersection.  A 
second  element  of  delay  (i.e.,  queueing  delay)  must  be  added  to  the  basic 
delay  (W)  caused  by  the  gating  of  vehicles.  This  delay  is  based  on  the 
average  queue  length  (Q) ,  which  is  given  by  the  equation: 

Q  =  V/(2*G*S/(R  +  G)  -  V) 

Then  the  average  delay  per  vehicle  is  the  delay  due  to  the  gating  effect 
of  the  signal  plus  the  time  spent  waiting  in  the  queue  of  length  Q,  which 
is  being  serviced  at  a  rate  (S  *  G)/(R  +  G) .  Combining  the  delay  per 
queued  vehicle  (which  is  the  reciprocal  of  the  service  rate)  with  the 
queue  length  and  combining  this  result  with  the  average  delay,  the  average 
delay  per  vehicle  including  delay  due  to  queueing  (D)  is : 

D  =  S*R*R/(2*(R  +  G)*(S  -  V))  +  V*(R  +  G) **2/ (2*S*G* (S*G  -  V*(R  +  G))) 

Tlie  computing  of  delay  at  an  intersection  requires  predetermination 
of  intersection  approach  capacity,  signal  phase  lengths,  and  cycle  times. 
The  BATS  Intersection  subroutine  recomputes  these  variables  for  each 
simulated  time  period  (usually  1  hour)  as  follows: 

•  Adjustments  to  intersection  approach  capacities  are  automatically 
made  based  on  the  number  of  left-turning  vehicles  at  an  inter¬ 
section  . 

0  Red  and  green  phase  lengths  as  well  as  cycle  times  are  auto¬ 
matically  computed  at  actuated  signals  [see  Webster  (1958)]. 

•  Left-turn  phases  may  be  specified  at  a  signalized  intersection, 
and  capacity  for  such  phases  may  be  input  or  automatically 
calculated  based  on  the  intersection  approach  capacity  for  the 
through  and  right-turning  traffic. 

Figure  10  is  a  flowchart  of  the  signalized  intersection  model. 


2. 5.2.2  Unsignalized  Intersections 

The  intersection  submodel  will  also  predict  delay  and  queue  length 
at  unsignalized  intersections.  As  with  a  signalized  intersection,  a 
queue  may  be  formed  because  of  random  arrivals  at  the  intersection. 
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FIGURE  10.  SIGNALIZED  INTERSECTION  MODEL  FLOW  DIAGRAM 


^  For  Each  Approach  J 

i 

>  For  Each  Turning  Movement  K 


J2  -  INDEX  OF  PHASE  GREEN  FOR  THIS  APPROACH 
AND  TURNING  MOVEMENT 


CAPlJ.K)  ■=  CAPlJ.KI  •  (PTIJ2I  -  CIII)/2l/COT 


VOLIJ.K)  *  1/COT  r  PT(J2I 
CAP  (J.KI 


NEWELL  MODEL  OF  DELAY  AND  QUEUE  LENGTH 

Wl  -  (COT  -  PT(J2)l^/(2  •  COT  *  (1  -  ylJ.KII) 

00  «  V0LIJ.KI/(2  •  (CAPlJ.KI  -  VOLIJ.KIII 
W3  =  DELAY  DUE  TO  QUEUE  FROM  PREVIOUS 
TIME  PERIOD 


W2  -  00  •  3600/CAP(J.K) 


DELA(L.KI  -  Wl  +  W2  ♦  W3 

WHERE  L  IS  THE 

LINK  ENTERING 

QQ(J,K)  ‘  00 

FROM  APPROACH  J 

MODIFIED  NEWELL  MODEL 
OF  DELAY  AND  QUEUE  LENGTH 
Wl  •  (COT  -  PT(J2ll/2 
OO  -  fIVOL.CAP.PTI 
Ql  -  VOL(J,K>  -  CAP(J,K)  ♦  QUEII.J) 
W3  -  TP 


QUEd.JI  >  QQ(J,U  r  QQ(J,3) 


If  V  Is  the  arrival  rate  of  vehicles  at  the  Intersection  and  C  is  capacity 
flow  rate  of  vehicles  departing  the  iiiierseit  Ion,  tin-  average  queue  lenp.th 
is: 


Q  =  v/(C  -  V) 

The  average  delay  experienced  by  vehicles  at  an  unsignalized  intersection 
is : 


D  =  V/(C*(C  -  V)) 

Figure  11  is  a  flowchart  of  the  unsignalized  intersection  model. 


2.5.3  Parking  Area 

The  EPA  Guidelines  (1978)  recommend  that  a  constant  be  added  to  the 
running  times  of  each  vehicle  operating  in  a  parking  lot  that  is  more  than 
80  percent  filled.  BATS  is  written  to  predict  air  pollution  in  a  more 
realistic  way.  First,  travel  time  is  increased  as  a  function  of  the 
percent  of  the  parking  lot  that  is  filled.  Second,  when  a  zone  is  filled 
to  capacity,  vehicles  that  would  have  traveled  to  this  zone  are  routed  to 
the  adjacent  zone  or  to  the  zone  with  the  next  highest  potential  for 
attracting  trips.  In  this  way,  parking  areas  are  filled  in  the  order  of 
their  potential  and  when  filled  operate  under  capacity-flow  conditions. 

The  parking  area  subroutine  predicts  running  time  for  the  vehicles 
arriving  or  leaving  each  parking  zone.  Running  time  is  divided  into 
three  parts: 

•  Time  to  travel  to  an  open  stall  starting  from  the  parking  zone 
entrance  or  time  to  travel  from  a  stall  to  the  exit  of  the  park¬ 
ing  zone 

•  Queueing  delay  caused  by  one  or  more  vehicles  waiting  for  a 
vehicle  to  back  out  of  a  stall 

•  Delay  to  vehicles  that  drive  around  a  parking  zone  looking  for 
an  empty  stall  either  because  of  an  unwillingness  to  park  in  a 
different  zone  or  because  there  are  no  other  alternative  zones  in 
which  to  park. 

A  flowchart  of  the  parking  lot  subroutine  is  shown  in  Figure  12.  Equa¬ 
tions  for  each  of  the  three  elements  of  running  time  are  as  follows: 

TT  ■=  PL/PLS 


where 

TT  is  travel  time  (sec)  to  vehicles  either  entering  or  leaving  the 
parking  zone 

PL  is  the  parking  lot  length  (ft) 

PLS  is  the  average  speed  (ft/sec)  of  vehicles  in  the  lot. 
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FtGURE  11.  UNSIGNAUZED  INTERSECTION  MODEL  FLOW  DIAGRAM 


DQ  =  MAX(PLBO**2/(2*(TP/VZ  -  Z))  ,  PLBO/Z) 

where 

VZ  is  the  number  of  arrivals  during  a  time  period 

DQ  is  the  delay  (sec)  caused  by  vehicles  backing  out  of  stalls 
during  near  capacity  conditions 

PLBO  is  the  time  (sec)  required  for  a  vehicle  to  back  out  of  and 
move  forward  enough  to  clear  a  parking  stall 

TP  is  the  length  of  a  time  period  (sec) 

AWQ  =  (PV  +  (VZ  -  ZV)/2  -  PC)*TP/ZV  [if  PV  ^  PC  and  PV  +  (VZ  -  ZV)/2  ^  PC] 

AWQ  =  (PV  +  VZ  -  ZV  -  PC)**2/(2*(VZ  -  ZV))*TP/ZV  [if  PV  +  (VZ  -  ZV)/Z  >  PC 
and  PV  <  PC] 

AWQ  =  (PV  -  PC)**2/(2*(ZV  -  VZ))*TP/ZV  [if  PV  >  PC  and  PV  +  VZ  -  ZV  <  PC] 

The  travel  time  equation  (TT  =  PL/PLS)  uses  PL  (parking  lot  length)  to 
represent  the  distance  a  vehicle  will  drive  when  it  is  traveling  through 
the  entire  parking  area.  PL  is  computed  by  finding  the  longest  side  of 
the  parking  zone,  assuming  the  parking  zone  is  rectangular,  finding  the 
shorter  side  of  the  rectangle,  and  adding  this  to  the  longer  side.  For 
example,  a  40,000-square  foot  parking  zone  with  a  longest  side  of 
200  feet  has  an  assumed  shorter  side  of  200  feet  and  a  PL  of  400  feet, 
while  a  zone  of  the  same  area  with  a  longest  side  of  400  feet  would 
have  a  100-foot  shorter  side  and  a  PL  of  500  feet. 

Vehicles  that  are  routed  to  each  parking  zone  may  be  of  six  types: 
autos,  light  duty  trucks  (<6000  pounds),  light  duty  trucks  (6,000-8,500 
pounds),  heavy  duty  trucks,  heavy  duty  diesel  (include  buses),  and 
motorcycles.  Autos  and  light  duty  trucks  routed  to  a  parking  area  or 
zone  will  enter  that  zone  and  incur  the  running  time  related  to  parking. 
Buses,  trucks,  and  other  vehicles  routed  to  a  parking  zone  will  not  look 
for  parking  places  nor  use  parking  stalls  in  the  parking  area. 
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SECTION  III 


DOCUMENTATION 


This  section  documents  each  subroutine  for  the  programmer  or  analyst 
who  is  using  the  BATS  program.  The  documentation  thus  is  useful  to  the 
computer  professional  who  wishes  to  change  a  subroutine  or  to  understand 
exactly  what  a  subroutine  does.  The  documentation  of  each  subroutine 
includes  a  description  of  the  purpose  of  the  routine,  a  data  description 
of  the  principal  parameters  used  by  each  subroutine,  the  reports  generated 
by  the  subroutine,  and  the  subroutines  called  by  this  subroutine.  The 
documentation  of  the  main  BATS  program  also  includes  a  flow  chart  showing 
the  major  subroutine  calls  made  by  the  BATS  model. 


3.1  BATS 


3.1.1  Purpose 

BATS  calls  the  major  functional  subroutines  and  performs  sequential 
iterations  of  the  modeling  process.  The  named  and  blank  common  areas 
are  all  defined  in  BATS  and  the  default  data  are  initialized,  A  flow¬ 
chart  of  the  BATS  program  is  given  as  Figure  13.  When  the  model  was 
debugged,  many  supplementary  FORTRAN  statements  were  included  in  the 
code  to  generate  additional  information,  and  these  statements  have  been 
made  into  comments  by  placing  an  in  column  1  on  the  printout.  The 
normal  comments  in  the  program  are  made  with  a  C  in  column  1.  A  complete 
program  listing  is  presented  in  Appendix  A. 


3.1.2  Data  Description 


Common  Name 

Description  of  Data 

/COMM/ 

Stores  the  general  information  about  the  run;  provides 
a  means  of  passing  parameters,  and  stores  the  itera¬ 
tion  variables  which  may  be  referenced  by  several 
different  subroutines. 

/CHD/ 

Holds  the  data  which  art'  output  as  the  first  line  of 
each  report. 

/LINK/ 

Stores  the  parameters  describing  each  link  of  the 
network;  data  are  generally  read  from  Card  Type  2. 

/INTRST/ 

Stores  the  parameters  describing  each  intersection  of 
the  network;  data  are  read  from  Card  Type  3. 

/ZONES/ 

Stores  the  parameters  describing  each  zone.  Some 
data  are  read  from  Card  Type  A;  other  values  such  as 
ZATTR,  ZGENR,  and  FCS  are  generated  by  subroutines. 

1 


INITIALIZE  TABLES  TO  ZERO  OH  DEFAULT  INPUT  VALUES 


INPUT  -  LINK,  INTERSECTION,  ZONE,  DEMOGRAPHIC 
VARIABLES,  TRIP  PURPOSE,  TRIP  GENERATION 
COEFFICIENTS,  GATE  COUNTS,  THROUGH  TRAFFIC 
COUNTS,  VEHICLE  TYPE  AND  LOAD  FACTORS 


DETERMINE  PARKING 
AREA  AND  TRIP 
DISTANCE 


END  OF  FILE  FOR> 
AQAM  OUTPUT 
-  TAPE  7  > 


CLOSE  PLOT  file 
ON  TAPE  27 


END 


3 


TRIPGEN  -  GENERATE  TRIP  PRODUCTIONS 
AND  ATTRACTIONS  IN  PERSON  TRIPS 
FOR  EACH  ZONE,  BY  TRIP  PURPOSE 


MNPATH  -  GENERATE  THE  MINIMUM  TRAVEL 
TIME  PATH  BETWEEN  ALL  GATES  AND  ZONES 


GRAVO  -  DISTRIBUTE  TRIPS  BETWEEN  ORIGIN 
AND  GATES,  GATES  AND  DESTINATIONS,  AND 
INTERNAL  ORIGIN  AND  DESTINATION  ZONES 


MODAL  -  USE  VEHICLE  TYPE  AND  LOAD  FACTOR 
TO  DETERMINE  VEHICLE  TRIPS  BETWEEN  ZONE-GATE 
PAIRS  AND  ZONE-ZONE  PAIRS.  DETERMINE 
BUS  LOAD  FACTOR  BASED  ON  EXPECTED  COST 
OF  TRAVEL  BY  BUS  OR  ALTERNATE  MEANS 


I 


FIGURE  13.  BATS  FLOW  CHART  SHOWING  MAJOR  SUBROUTINES  CALLED 
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rSMOOTH  -  THE  DATA  IS  STORED  AS  ZONE-GATE  VEHICLE  VOLUMES  FOR' 
EACH  GATE  AND  ZONE.  THESE  ARRAYS  ARE  FACTORED  TOWARD  THE 
GATE  COUNTS,  A  USER  OPTION  ALLOWS  THE  GATE  COUNTS  TO  REMAIN 
CONSTANT  AND  THE  ZONE  COUNTS  FACTORED  (BASE  CASE)  OR  THE 
ZONE  COUNTS  CONSTANT  AND  THE  GATE  COUNTS  FACTORED 
(PREDICTIVE  CASE).  REDISTRIBUTE  VEHICLE  TRIPS  FROM 
V OVER  CAPACITY  PARKING  ZONES  TO  UNDER-CAPACITY  PARKING  ZONES  J 


ASSIGN  -  ASSIGN  MILITARY  AND  CIVILIAN  TRIP  PURPOSE 
VEHICLE  TRIPS  TO  THE  NETWORK.  MAINTAIN  VEHICLE 
TYPE  AND  HOT/COLD  START  COUNTS  ON  EACH  LINK 


MODEL  LINK  SPEEDS,  INTERSECTION  DELAYS, 
PARKING  LOT  TRAVEL  TIMES 


PRINT  MODEL  RESULTS 


AQAM  OUTPUT 
REQUIRED 


AQAM  -  SAVES  A  FILE  THAT  WILL  BE  READ  BY  AQAM. 

THIS  HAS  SPEED,  VEHICLE  TYPE,  VOLUME  HOT/COLD  START  AND 
MILITARY  VEHICLE  SPEED,  AND  TYPE,  VOLUME  DATA  FOR  EACH 
LINE  SOURCE.  THE  FILE  IS  IN  UNITS  OF  YEARLY  VOLUME, 
FACTORS  ARE  THE  INVERSE  OF  THOSE  FOUND  IN  AQAM  FOR 
CORRECTING  TO  HOURLY  VALUES  FROM  YEARLY  VALUES 


PLOT  OUTPUT 
REQUIRED 


PLOT  -  OUTPUTS  A  PLOT  FILE  INCLUDING  x,  y 
COORDINATES  OF  EACH  LINK,  VOLUME,  CAPACITY 
ZONE  DEFINITIONS  (SELECT  TIME  PERIOD  TO  BE  PLOTTED) 


FIGURE  13.  BATS  FLOW  CHART  SHOWING  IVIAJOR  SUBROUTINES  CALLED  (Concluded) 
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/DEMVAR/ 


Stores  the  parameters  to  be  used  to  generate  trip 
ends  for  each  zone;  data  are  read  from  Card  Types  5 
and  15. 


/TRIP/ 

/GATE/ 

/VOLUME/ 

/VEHLOD/ 

/TROUT/ 

/SHIFT/ 

/PARKZ/ 

/RESET/ 

/CAPMAX/ 


Stores  the  coefficients  and  parameters  used  to  generate 
trips  from  demographic  variables;  data  are  read  from 
Card  Type  6.  The  data  are  used  principally  by  sub¬ 
routine  TRIPGEN. 

Stores  the  gate  counts  and  other  gate  definition 
data  read  from  Card  Type  7. 

Stores  the  through  traffic  counts  read  from  Card 
Type  8. 

Stores  the  data  on  vehicle  load  factors  for  each  type 
of  vehicle  as  read  from  Card  Type  9. 

Stores  the  truck  route  data  as  read  from  Card  Type  10. 

Stores  the  shift  count  numbers  as  read  from  Card 
Type  11. 

Stores  the  data  used  to  determine  parking  lot  running 
times  by  subroutine  PARKING;  data  are  generated  from 
elements  of  LINK,  ZONE,  and  TRIP  conmon  blocks. 

Stores  the  results  of  the  simulation;  data  are  used  to 
initialize  the  program  on  subsequent  iterations. 

The  links  with  volumes  that  exceed  a  capacity  will 
be  specially  marked  on  the  plot  output. 


3.1.3  Reports 

I.l.  NETWORK  SUMMARY  PARAMETERS  FOR  TIME  PERIOD.  This  totals 

delay,  vehicle  miles  and  other  parameters  on  all  links  of  the 
network.  This  report  gives  a  general  measure  of  the  effective¬ 
ness  of  the  traffic  network. 


3.1.4  Subroutines  Called 

The  following  subroutines  are  called:  INPT,  ZAREA,  MNPATH,  TRIPGEN, 
GRAVO,  MODAL,  SMOOTH,  ASSIGN,  INSEC,  INSECU,  COORXY,  PARKING,  PLOTI, 
PLOTA,  PLOTP,  AQAMQ.  Other  subroutines  (described  in  Sections  3.18-3.37) 
are  called  by  these  major  subroutines. 


3.2  INPT 
3.2.1  Purpose 

This  routine  reads  all  data  on  Card  Types  0-16  that  are  used  to  run 
an  Iteration  of  the  BATS  model.  Some  manipulation  of  the  Input  data  is 
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undertaken  in  INPT;  for  instance,  x,y  coordinates  of  link  end  points 
are  converted  to  UTM  (Universal  Transverse  Mercator)*  coordinates.  All 
input  cards  are  printed  out  after  being  read,  to  tell  the  user  where  an 
input  error  occurred  and  to  permanently  record  the  input  data  used  to 
produce  the  resulting  reports. 


3.2.2  Data  Description 

Table  1  describes  all  input  variables,  by  input  card  type,  field, 
format,  symbol,  dimensional  units,  value  limits,  default  values,  and 
variable  description.  Instructions  for  assembling  the  input  data  are 
contained  in  the  "BATS — Field  Data  Collection  and  Reduction  Guide" 
(Swope,  et  al.,  1979). 


3.2.3  Reports 

The  input  data  are  listed  using  the  same  format  specifications  as 
used  to  read  input  cards. 


3.2.4  Subroutines  Called 

The  following  subroutines  are  called:  AQAMl,  COORXY,  LETTER,  PLOT, 
PLOT!. 

3.3  ZAREA 


3.3.1  Purpose 

This  routine  will  compute  the  area  of  a  3  to  12  sided  zone.  It  uses 
the  first  side  of  the  zone  to  construct  a  trapezoid  using  the  X  axis  and 
two  sides  parallel  to  the  Y  axis.  Then  it  uses  the  second  side  to  con¬ 
struct  another  trapezoid,  and  so  forth,  adding  or  subtracting  the  areas 
of  the  trapezoids  created.  This  routine  was  part  of  the  software  on  a 
programmable  calculator,  and  is  explained  here  more  completely: 

(1)  The  trapezoid  rule  states  the  area  of  a  trapezoid  is  the  total 
length  of  the  parallel  sides  times  the  distance  between  them 
divided  by  2 . 

(2)  If  and  X2,Y2  are  the  end  points  of  a  line  segment  used 

to  define  an  area,  then  the  area  of  the  trapezoid  with 


See  "AQAM — Field  Data  Collection  Guide"  (1975),  page  10,  for  information 
on  finding  UTM  coordinates  of  a  point. 
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Table  1  (Continued) 
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Table  1  (Continued) 


BATS  BASIC  INPUT  INFORMATION  CARDS  0  THROUGH  16 


Table  1  (Continued) 
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(X],Yj;  X2,Y2).  (Xi.Yj;  Xj^.O),  (X2,Y2;  X2,0),  and 
(X|,0;  X2»0)  a«  tht  four  sldc-.s  is: 

(Y^  +  Y^XX^  - 
(3)  The  area  of  the  zone  is  then: 

A  =  ((Y.  +  Y„)(X,  -  X„)  +  (Y„  +  Y„)(X„  -  X.)  +  ...  (Y  +  Y  ) (X  ,  -  X  ) 
\1  zi  t  L  i  L  i  n-i  nn-i  n 

+  (Y  +  Y  )(X  -  X  ))/2 

n  In  If 

Note  the  X]^,Yl»  X2,Y2»  ^3X3*  •••  Xn>Yn  n>ust  be  sequential 
points  that  define  the  perimeter  of  the  zone. 


3.3.2  Data  Description 


Xi(L),Yi(L),X2a),Y2(L) 

PNOS(Z) 

ZLINKSd.Z) 

PX(NL),PY(NL) 


TRIARE(NL) 


Meaning 

The  x,y  coordinates  of  link  L 

The  number  of  links  defining  zone  Z. 

The  link  numbers  (1=1,12)  defining 
the  area  Z. 

The  coordinates  of  the  NLth  point 
defining  a  zone. 

Tlie  number  of  points  defining  a  zone. 
The  area  of  the  NLth  trapezoid. 


3.3.3  Reports 

When  the  x,y  coordinates  of  a  link  do  not  indicate  connection  to 
the  previous  link,  the  error  message,  "ZONE  Z  is  disconnected  between 
links  L  and  LI,"  is  printed  out.  The  area  of  this  zone  will  be  incor¬ 
rectly  computed  based  on  the  sides  that  were  connected. 


B.l.  ZONE  PARKING  CAPACITIES  AND  TRIP  LENGTHS  is  a  table  of  the 
results  of  ZAREA. 


3 . 4  MNPATH 


3.4.1  Purpose 

Tills  subroutine  is  responsible  for  generating  the  minimum  path 
between  each  origin  and  destination  zone  in  the  network.  The  method 
generates  the  minimum  path  from  one  origin  AA  to  all  links  of  the  network. 
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(Recall  that  a  link  is  one  direction  of  flow  on  a  street  and  has  a  load 
point  halfway  between  each  end.)  Two  arrays  are  initialized  before 
calling  MNPATH.  The  TT(L,J)  array  represents  the  travel  time  from  link 
L  to  a  connected  link  J.  The  C(L)  array  represents  the  cost  of  travel 
from  the  origin  zone  to  any  link  L  of  the  network.  The  C(L)  array  is 
initialized  to  zero  for  any  link  connected  to  the  origin  zone  and  to 
infinity  for  any  iinronn(“c  ted  lirjk. 

Tile  algorillim  described  below  is  carried  out  using  these  two  arrays. 
On  completion,  the  (1(1.)  array  represents  the  minimum  travel  time  (or 
cost)  from  the  origin  zone  to  each  link  I,  of  tlu^  network.  During  this 
process,  the  K(l,)  array  is  generated,  which  contains  the  Jink  number 
th.it  immed  i.itel  y  precedes  link  1.  in  the  route  from  the  origin  to  L. 


A  Minimum  Path  Algorithm  for  Route  Generation: 

Step  1 — Initialize  TT(L,J)  with  the  constant  cost  in  traveling 
from  link  L  to  three  possible  adjacent  links,  J  =  J2,  ^3. 

Step  2 — For  an  origin  zone  initialize  with  a  large  number  the 
cost,  C(L),  to  travel  from  the  origin  to  each  link  in  the 
network  (C(L)'  “  for  all  L) . 

Step  3 — For  the  links  adjacent  to  the  origin  zone  store  a 
travel  time  ca)st  In  G(l.)  (I.  =  1,^,  h2»  •••> 

Step  /\ — Set  the-  travel  time  costs  of  Step  3  negative,  to  flag 
th.iL  .ill  links  connected  to  these  may  have  an  improved  travel 
Lime  cost. 

Step  5 — Kind  a  negative  C(L)  and  the  connected  links  J2,  J3 
to  link  L.  If  no  negative  C(L)  exists  then  do  Step  9.  If  no 
connected  links  exist  then  do  Step  8. 

Step  6 — For  each  J  =  J]^,  J2,  J3  test  |C(L)1  +  TT(L,J)  <  |C(J)|. 
If  the  relationship  does  not  hold  true  then  do  Step  7.  Other¬ 
wise  replace  C(J)  by  |C(L)|  +  TT(L,J)  and  set  C(J)  negative. 

The  route  is  saved  by  storing  the  previous  link  to  J  (i.e., 

R(J)  =  L). 

Step  7 — Repeat  from  Step  6  until  all  J's  have  been  tested. 

Step  8 — Set  C(L)  to  jc(L)|  and  repeat  from  Step  5. 

Step  9 — The  C(L)  array  is  the  cost  of  travel  from  the  origin 
to  each  link  h.  Repetition  from  Step  2  will  generate  new 
costs  of  travel  from  other  origins. 
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3.4.2  Data  Description 


Symbol 

C(L) 

R(L) 

TT(L,J) 

LINKS (I) 
NL 

Cl 

C2 


Meaning 

Cost  of  time  going  from  origin  zone  to  link  L 

Link  previous  to  link  L  in  route  from  origin  zone  to 
link  L 

Travel  time  going  from  link  L  straight  (J=l),  right  (J=2) , 
or  left  (J=3) 

The  links  accessing  the  zone  (or  gate)  of  origin 

The  number  of  links  accessing  the  zone  (or  gate)  of 
origin 

The  cost  to  get  from  the  center  of  the  origin  zone  to 
LINKS(l) 

The  cost  to  get  from  the  center  of  the  origin  zone  to 
LINKS(2),  LINKS(3),  etc. 


3.5  TRIPGEN 
3.5.1  Purpose 

The  trip  generation  routine  generates  trip-end  productions  and 
attractions  as  linear  functions  of  up  to  10  demographic  variables  for  up 
to  15  trip  purposes.  The  coefficients  used  to  predict  linear  productions 
and  attractions  are  specified  by  the  user,  or  they  may  be  generated  by 
default  using  a  resident  array  of  trips  related  to  land  use. 


3.5.2  Data  Description 


Symbol  Meaning 

Initial  (COMMON/TRIP/) 

VAR(J,Z)  The  Jth  demographic  variable  associated  vith 

zone  Z. 


COEFO(I,K) 

COEFDd.K) 

NVAROd.K) 

NVARD(I,K) 

LANDOd.K) 


The  Ith  coefficient  for  generating  productions 
for  trip  purpose  K. 

The  Ith  coefficient  for  generating  attractions 
for  trip  purpose  K. 

The  demographic  variable  index  associating 
the  Ith  production  coefficient  with  a  demo¬ 
graphic  variable. 

The  demographic  variable  index  associating 
the  Ith  attraction  coefficient  with  a 
demographic  variable. 

The  land  uses  associated  with  trip  purpose 

K. 
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Initial  (COMMON/ TRIP/)  (continued) 


PLUALU(L,LU,I) 

The  number  of  trips  from  land  use  L  zones  to 
land  use  LU  zones  in  the  a.m.  off-peak, 
noon  rush,  or  p.m.  off-peak  hours  (1=1,2,  or 

3,  respectively) . 

NTRIP 

The  number  of  trip  purposes. 

NTRIPC 

The  number  of  trip  purposes  using  civilian 
vehicles. 

LANDU(Z) 

The  land  use  designated  for  zone  Z. 

KCSP(K,1) 

Fraction  of  cold  start  vehicxes  with  external 
origins  for  purpose  K. 

FCSP(K,2) 

Fraction  of  cold  start  vehicles  with  internal 
origins  for  purpose  K. 

Intermediate 

ALUA(L,LU) 

The  number  of  trips  going  from  land  use  L 
zones  to  land  use  LU  zones  in  the  time  period 
being  simulated. 

CSPLU(LU) 

The  column  sum  (attractions)  of  the  PLUALU 
array  for  zones  with  land  use  LU. 

RSPLU(LU) 

The  row  sum  (productions)  of  the  PLUALU 
array  for  zones  with  land  use  LU. 

CSUM(LU) 

The  column  sum  (attractions)  of  ALUA  array 
for  zones  with  land  use  LU. 

RSUM(LU) 

The  row  sum  (productions)  of  the  ALUA  array 
for  zones  with  land  use  LU. 

VARZDV 

The  weighted  average  of  all  demographic 
variables  associated  with  a  zone. 

SFTO 

Sum  of  the  flex  time  trips  from  all  origins. 

SFTD 

Sum  of  the  flex  time  trips  to  all  destinations. 

TSHO 

Total  shift  trips  from  all  origins. 

TSHD 

Total  shift  trips  to  all  destinations. 

TGATEC 

Total  gate  count  for  four  15-minute  periods. 

GAFRO 

Fraction  of  trips  arriving  in  15  minutes 
from  all  origins. 

GAFRD 

Fraction  of  trips  arriving  in  15  minutes  to 
all  destinations. 

MATRIX(G,Z) 

Matrix  associating  gate  G  with  zone  Z 
(1  =  associated;  0  =  not  associated). 

FRAC 

Fraction  of  all  flex  pool  trips  to  be 
routed  to  a  zone. 
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Intermediate  (Continued) 

SHFPCT(Z,1,T)  Shift  employee  origins  from  zone  Z,  time 

period  T. 

SHFPCT(Z, 2, T)  Shift  employee  destinations  to  zone  Z, 

time  period  T. 


The  number  of  person  trips  originating  at 
zone  Z  for  trip  purpose  K. 

The  number  of  person  trips  destined  for 
zone  Z  for  trip  purpose  K. 

Fraction  of  cold  starts  for  zone  Z. 


Resultant 

NTO(Z,K) 

NTD(Z,K) 

FCS(Z) 


3.6  GRAVO 


The  GRAVO  subroutine  distributes  the  person  trips  associated  with 
each  zone  into  person  trips  associated  witli  each  pair  of  zones  and  a 
gate.  A  gravity  model  is  used  to  distribute  trips  among  zone  pairs  based 
on  zone-gate  and  gate-zone  impedances  (usually  travel  time)  and  zonal 
mass  (number  of  trip  ends).  The  general  form  of  the  model  is  as  follows 
where  Z  is  an  external  zone  and  Z1  is  an  interal  zone: 
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OV(Z,G)  =  VD(Z1,G)  +  GCNT15(G,2, I)/(CGO(G)  +  CGD(G)**3) 

all  Z1 

*NT0(Z,K)*NTD(Z1,K)/  NTD(Z1,K) 


The  equation  makes  the  number  of  persons  going  from  zone  Z  to  zone  Z1 
through  gale  G  proportional  to  the  gate  count  at  gate  G  and  the  number 
of  person  trips  being  made  between  zone  Z  and  Z1  Indirectly  proportional 
to  the  travel  time  to  the  gate,  plus  the  cube  of  the  travel  time  within 
the  base. 


3.6.2  Data  Description 


Symbol 

Initial 

Meaning 

NTO(Z,K) 

The  number  of  person  trips  originating  at  zone  Z 
for  trip  purpose  K 

NTD(Z,K) 

The  number  of  person  trips  destinated  to  zone  Z 
for  trip  purpose  K 

GCNT15(G,1,I) 

The  number  of  vehicles  exiting  through  gate  G 
during  time  period  I 

GCNT15(G,2,I) 

The  number  of  vehicles  entering  through  gate  G 
during  time  period  I 

CG(G, L) 

The  cost  to  travel  from  gate  G  to  link  L  (also 
assumed  to  be  the  cost  of  travel  from  link  L  to 
gate  G) 

Intermediate 

SGCOl 

The  sum  of  the  existing  gate  counts  for  this  time 
period 

SGC02 

The  sum  of  the  entering  gate  counts  for  this  time 
period 

TOSUM(K) 

The  total  trips  originating  on  base  for  purpose  K 

TDSUM(K) 

The  total  trips  destined  on  base  for  purpose  K 

ZGT 

The  trips  originating  at  an  internal  zone  Z  for 
all  trip  purposes 

ZAT 

The  trips  destined  for  an  internal  zone  Z  for  all 
trip  purposes 

Resultant 

OV(Z,G) 

The  number  of  person  trips  going  from  zone  Z  to 
gate  G. 

VD(Z,G) 

The  number  of  person  trips  going  to  zone  Z  from 
gate  G. 
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Resultant  (Continued) 


NTO(Z,K) 

NTD(Z,K) 


Tie  number  of  person  trips  going  from  internal 
zone  Z  to  another  internal  zone  for  purpose  K. 

The  number  of  person  trips  going  to  internal 
zone  Z  from  another  internal  zone  for  purpose  K. 


3.6.3  Reports 

D.l.  ORIGIN  TO  GATE  AND 
D.2.  ORIGIN-DESTINATION 
D.3.  ORIGIN-DESTINATION 
D.4.  ORIGIN-DESTINATION 

3.7  MODAL 


GATE  TO  DESTINATION  TRIPS 
ARRAY 

ARRAY  FOR  CIVILIAN  VEHICLE  TRIPS 
ARRAY  FOR  MILITARY  VEHICLE  TRIPS 


3.7.1  Purpose 

This  subroutine  determines  the  mode  of  travel  of  all  person  trips 
going  between  all  origin  and  destination  zones.  The  0V(1Z,IG)  and 
VD(IZ1,IG)  arrays  store  the  number  of  persons  originating  at  zone  IZ  and 
going  through  gate  IG,  and  the  number  of  persons  going  through  J(;  and 
destined  for  zone  IZl.  All  these  person  trips  must  be  assigned  to  some 
mode  of  travel,  which  is  one  of  six  types  of  civilian  vehicles  and 
seven  types  of  military  vehicles.  The  percentage  and  the  load  factor 
for  each  type  of  vehicle  for  each  zone  are  known.  Thus,  the  number 
of  vehicle  trips  between  IZ  and  IG  and  between  IG  and  IZl  is  typically 
easy  to  calculate. 

The  exception  is  military  vehicle  type  7 — a  bus.  A  bus  may  not 
have  a  load  factor,  in  which  case  it  must  be  computed  in  the  MODAL 
subroutine. 


3.7.2  Data  Description 


Symbol 

OV(Z,G) 

VD(Z,G) 

VEHTYP(1,Z) 

VTYPM(1,Z) 


Meaning 

Origin  zone  Z  to  gate  G,  demand  volume  in  persons. 

Gate  G  to  destination  zone  Z,  demand  volume  in  persons. 

Percent  of  civilian  vehicles  of  each  type  I,  at  each 
zone  Z. 

Percent  of  military  vehicles  of  each  type  1,  at  each 
zone  Z.  For  VTYPM(7,Z)  this  is  the  actual  number 
of  buses. 
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FRAMIL(Z,  I) 


VLF(I) 

VLFM(I) 

0V(Z,NG1) 

VD(Z,NG1) 

TP 

VZ(Z) 

ZV(Z) 

NEXT 

NZONES 

NGATE 


Fraction  of  person  trips  that  use  military  vehicles 
for  each  zone  Z;  FRAMIL(Z,1)  is  the  fraction  that 
originates  at  a  zone  Z,  and  FRAMIL(Z,2)  is  the  frac¬ 
tion  destined  for  a  zone  Z. 

Load  factor  in  persons/vehicle  for  civilian  vehicles 
of  type  I. 

Load  factor  in  persons/vehicle  for  military  vehicles 
of  type  I. 

Internal  origin  zone  Z  productions  of  internal  trips. 
Internal  origin  zone  Z  attractions  of  internal  trips. 
Length  of  time  of  simulation. 

Number  of  parking  trips  to  each  zone  Z. 

Number  of  parked  vehicles  leaving  each  zone  Z. 

Number  of  exterior  zones. 

Number  of  interior  and  exterior  zones. 

Number  of  gates  IG. 


'i .  7 .  '3  Report s 

E.l.  MODAL  SPLIT  VEHICLE  LOAD  FACTORS 

E.2.  ORIGIN  TO  GATE  AND  GATE  TO  DESTINATION  TRIPS 


3 . 8  SMOOTH 


3.8.1  Purpose 

The  gate  counts  can  be  used  to  calibrate  the  vehicle  trips  predicted 
between  zones.  The  external-to-internal  and  internal- to- external  trips 
should  equal  the  gate  counts.  Thus,  the  sura  of  external  trip  productions 
can  be  made  to  equal  the  sum  of  the  vehicles  coming  through  all  the  gates, 
and  the  internal  to  external  trip  productions  can  be  made  to  equal  the 
sum  of  the  outgoing  gate  counts.  Another  calibration  function  is  per¬ 
formed  to  route  vehicles  through  the  appropriate  gates.  So  far  the 
model  has  weighted  the  trips  predicted  through  each  gate  by  the  gate 
counts.  A  matrix  multiplication  is  undertaken  first  to  make  the  columns 
of  the  OV(Z,G)  and  VD(Z,G)  arrays  add  up  to  the  gate  counts  and  second 
to  make  the  rows  of  the  OV  and  VI)  arrays  add  up  to  zone  productions  and 
attracLio ns. 


3.8.2  Data  Desc. rij[)l  i_on 

Symbo 1  Mean  I n^ 

OV(Z,G)  Number  of  vehicle  trips  going  from  zone  Z  to  gate  G. 
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VD(z,(;) 

fiCSUMl 

OVSUMl 

FEXGEN 

FEXATT 

FINGEN 

FIN  ATT 

ZGENR(Z) 

ZATTR(Z) 

SOVIZ 

SVDIZ 

SOVIG 

SVDIG 

ZV(Z) 

VZ(Z) 

PLA(Z) 

PD 

PV(Z) 


Number  of  vehicle  trips  going  to  zone  Z  from  gate  G. 

Sum  of  gate  counts  entering  the  base. 

I'reclicted  count  of  entering  vehicles. 

The  factor  for  converting  GCSUMl  to  OVSUMl. 

The  factor  for  converting  predicted  exiting  vehicles 
to  exiting  gate  counts. 

The  factor  for  converting  predicted  exiting  vehicles 
to  gate  counts. 

The  factor  for  converting  predicted  entering  vehicles 
to  gate  counts. 

The  total  generations  from  zone  Z  that  go  through  a  gate. 

The  total  attractions  from  zone  Z  that  go  through  a  gate. 

Sum  of  the  OV  array  for  all  zones  and  one  gate. 

Sum  of  the  VD  array  for  all  zones  and  one  gate. 

Sum  of  the  OV  array  for  one  zone  and  all  gates. 

Sum  of  the  VD  array  for  one  zone  and  all  gates. 

The  number  of  vehicles  that  park  in  zone  Z. 

The  number  of  vehicles  that  leave  from  parking  places  in 
zone  Z. 

The  parking  area  of  a  zone. 

Tfie  parking  density  of  a  zone. 

The  number  of  parked  vehicles  in  a  zone. 


3.8.3  Reports 

F.l.  CALIBRATION  FACTORS 

F.2.  ORIGIN  TO  GATE  AND  GATE  TO  DESTINATION  TRIPS. 


3.9  ASSIGN 

3.9.1  PuVpose  — . . 

This  subroutine  assigns  vehicle  trips  as  stored  in  the  NTO,  NTD, 

OV,  and  VD  arrays  to  the  links  of  the  network,  and  maintains  vehicle 
counts  on  each  link  for  each  through,  right,  left,  and  terminating 
movement,  and  for  civilian  types  1-6,  military  types  1-6,  hot  start, 
and  cold  start  vehicles.  The  subroutine  depends  on  the  MNPATH  subroutine 
to  return  the  minimum  time  path  between  zones.  Two  paths  between  each 
gate-zone  pair  or  internal  zone-internal  zone  pair  are  found.  Trips  are 
allocated  to  the  alternate  routes  based  on  travel  time.  Vehicle  type  data 
are  maintained,  as  are  hot/cold  start  data. 


65 


3.9.2  Data  Description 


Symbol 
Initial 
OV(Z,G) 
VD(Z,G) 
NTO(Z,K) 
NTI)(Z,K) 
!<(;((;,  1,) 

CG((;,  I.) 

LCON(L,J) 

ZLINKS(I,Z) 

Intermediate 

ZVEH(I) 

ZVEHM(I) 

DIV 

ZTRIP 

ZTRIPl 

OD 

Resultant 
COUNT (1,L) 
COUNT ( 2, L) 
COUNT ( 3, L) 
COUNT ( 4, L) 
COUNT ( 5, L) 
COUNT (6, L) 
COUNT(7,L) 
COUNT(8,L) 
COUNT (9, L) 
COUNT ( 10, L) 
COUNT (11,L) 
COUNTf  12,  I,) 
GOUN'I'fl  3,  l.j 
COUNIT  14,1.; 
COUNT  (lb,  I,; 


Meaning 

Vehicles  going  from  zone  Z  to  gate  G 

Vehicles  going  from  gate  G  to  zone  Z 

Vehicles  going  from  zone  Z  to  internal  zones 

Vehicles  going  to  zone  Z  from  Internal  zones 

Route  from  gate  G  to  each  link  L 
Coat  from  gate  G  to  each  link  L 
Link  connecting  to  L  going  J  direction 
Links  bordering  zone  Z,  (1=1  to  12) 


Civilians  vehicles  type  I — used  by  SUMIT 
Military  vehicles  of  type  I — used  by  SUMIT 
Divisor  applied  to  trips  in  SUMIT 
Trips  from  zone  Z 
Trips  to  zone  Z1 

Trips  going  from  zone  Z  to  zone  Z1 


Vehicle  count  going  through  on  link  L 
Vehicle  count  going  right  on  link  L 
Vehicle  count  going  left  on  link  L 
Vehicle  count  terminating  on  link  L 


Civil Ian 

vehicles 

of 

type 

1 

.  Civilian 

vehicles 

of 

type 

2 

Civilian 

vehicles 

of 

type 

3 

Civilian 

vehicles 

of 

type 

4 

Civilian 

vehicles 

of 

type 

5 

Civilian 

vehicles 

of 

type 

6 

Military 

vehicles 

of 

type 

1 

Mil itary 

vehicles 

of 

type 

2 

Mill Lnry 

veil  Icl  es 

of 

type 

3 

Mill  La  ry 

viTi  1  clf.s 

(>l 

type 

4 

Military 

vc'h  1  cles 

of 

type 

5 

hf, 


Resultant  (Continued) 

C0UNT(16,L)  Military  vehicles  of  type  6 
C0UNT(17,L)  Nunber  of  cold  start  velilcles 
COUNT ( 18, L)  Number  of  hot  start  vehicles 

3.9.3  Reports 

G.l.  ASSIGNMENT  COUNTS  AND  ASSOCIATED  COMPUTER  RUN  TIME 
G.2.  VEHICLE  COUNT,  TYPE,  AND  HOT/COLD  STARTS 

3.9.4  Subroutines  Called 

The  following  subroutines  are  called:  MNPATH,  SUMIT. 


3.10  INSEC 
3 . 10 . 1  Purpose 

This  subroutine  computes  delays  and  queue  lengths  at  signalized 
intersections  within  the  simulated  network.  Three  types  of  signal  con¬ 
trollers  can  be  specified.  The  first  type  is  a  fixed-time  controller  in 
which  the  north-south,  east-west,  and  left-turning  phases  are  all  of  fixed 
time  length.  The  other  types  of  controllers  are  vehicle-actuated,  tor 
which  types  the  INSEC  routine  must  determine  the  length  of  each  signal 
phase.  Phase  lengths  are  determined  based  on  volume  to  green  capacities 
of  each  approach  to  an  intersection  (Webster,  1958). 

Delays  and  queue  lengths  are  determined  at  an  intersection  on  the 
basis  of  phase  time,  approach  volume,  and  capacity  (Newell,  1965).  The 
referenced  method  of  determining  delay  assumes  a  period  of  time  during 
which  traffic  demand  remains  constant,  and  does  not  Include  a  means  of 
estimating  queueing  when  volume  exceeds  capacity.  For  BATS,  a  method 
was  devised  to  provide  continuous  functions  for  delay  and  queue  length 
when  volumes  approach  or  exceed  the  intersection  approach  capacity. 

This  method  establishes  a  constant  rate  of  increase  in  queue  length  beyond 
the  point  where  average  delay  at  the  intersection  equals  one  cycle  length. 
The  queue  is  assumed  to  Increase  at  a  constant  rate  beyond  this  point 
until  volume  exceeds  capacity,  at  which  point  the  queue  increases  at  the 
rate  that  volume  exceeds  capacity.  A  queue  is  dissipated  at  the  rate  that 
capacity  exceeds  volume  until  the  average  queue  length,  based  on  the 
volume  to  capacity  ratio  of  the  intersection  approach,  is  reached. 

The  results  of  the  model,  stored  in  the  DELA  and  QUE  arrays,  are 
delay  and  queue  length  for  each  approach  to  an  intersection. 
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'i.10.2  IJata  Ucisc rXptXon 


Symbol 

N1 

LIN(I,4) 

ITYPC(I) 

ICYCL(I) 

PH(I,4) 

CI(I) 

GCAP(I,4) 

QUE(I,4) 


Meaning 

Number  of  intersections 

Links  approaching  intersection  I  from  the  north,  east, 
south,  and  west 

Type  of  control  for  intersection  I 
Cycle  time  of  intersection  I 

Phase  time  on  the  north-south,  east-west,  north-south 
left  turn,  or  east-west  left  turn  to  intersection  I 

Clearance  interval  (i.e.,  yellow  duration)  at  intersec¬ 
tion  I 

Green  capacity  on  the  north,  east,  south  and  west  approach 
to  intersection  I 

The  average  queue  length  at  intersection  I 


3.10.3  Reports 

H.2.  INTERSECTION  DELAYS  AND  QUEUEING 


3.11  INSECU 


3.11.1  Purpose 

This  subroutine  predicts  the  delay  and  queue  length  at  unsignallzed 
intersections  in  the  network.  Three  types  of  intersections  can  be 
handled:  uncontrolled,  two-way  stop,  and  four-way  stop.  If  not 
specified  as  an  input  parameter,  the  INSECU  routine  determines  an  inter¬ 
section  capacity  based  on  the  Highway  Capacity  Manual  (HCM)  (1965). 

For  two-way  stops  and  uncontrolled  intersections,  the  HCM  recommends 
that  capacity  be  computed  as  if  a  signal  were  present  and  as  if  the  signal 
split  ratio  equaled 


VolumCj^  Width^ 
Volume^  Width^ 


At  four-way  stops,  capacity  is  a  function  of  the  number  of  lanes  and 
the  demand  split  among  approaches,  as  shown  in  Tables  6.7  and  6.8  of 
the  Highway  Capacity  Manual.  Queue  length  at  an  intersection  is  computed 
from  classical  queueing  theory,  i.e.,  QO  =  l/(CAP/VOL)  -  1.  Average 
delay  is  a  function  of  queue  length  times  the  time  to  process  each 
vehicle. 
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3.11.2  Data  Description 


Symbol 

NSTOPS(L) 

PNS 

DELA(L, J) 
QUEd.J) 


Meaning 

Number  of  stops  made  on  I  ink  I. 

Proportion  of  stops  on  a  link 

Delay  going  through,  right,  and  left  on  Jink  I. 

Queue  length  of  vehicles  at  intersection  I,  turning 
movement  J. 


3.11.3  Reports 

H.2.  INTERSECTION  DELAYS  AND  QUEUEING 


3.12  PARKING 
3.12.1  Purpose 

This  routine  models  parking  zone  flow  and  determines  vehicle 
running  time  for  a  zone.  The  running  time  of  an  average  vehicle  is 
considered  to  consist  of  three  elements:  (1)  the  time  to  travel  from 
the  edge  of  the  zone  to  a  stall  or  the  time  to  travel  from  a  stall  to 
the  edge  of  a  zone  (symbolized  by  TTA  and  TTD,  respectively);  (2)  the 
delay  experienced  by  arriving  vehicles  while  waiting  for  vehicles  to 
back  out  of  stalls  (DQ) ;  (3)  the  average  wait  in  a  queue  until  a  parked 
vehicle  leaves  the  zone  (AWQ). 

TTA  is  computed  on  the  basis  of  the  length  and  speed  of  an  average 
trip  into  the  parking  zone.  It  would  appear  that  the  length  and  speed 
of  an  average  trip  vary  in  direct  and  indirect  proportion,  respectively, 
to  the  utilization  of  a  zone,  and  preliminary  analysis  jf  available 
data  tends  to  support  this  supposition.  TTD  is  the  same  as  TTA,  but 
with  a  time  added  to  back  out  of  a  parking  stall  (PLBO) . 

DQ  is  the  queued  delay,  due  to  interrrupted  flow  lasting  PLBO 
seconds,  experienced  by  vehicles  arriving  in  the  parking  zone.  AWQ  is 
a  minimum  value  or  is  the  time  required  to  service  each  vehicle  times 
an  average  queue  length  estimated  from  the  utilization  of  the  parking 
zone. 


3.12.2  Data  Description 


Symbol  Meaning 

PC  Parking  lot  capacity 

PVZ  Parked  vehicles  at  end  of  time  period 

PLS  Parking  lot  speed 


AWQ 

TTA 

IXJ 

ZV(Z) 

VZ(Z) 


Average  wait  in  queue 
Travel  time  arriving 

Uacking  queue-waiting  for  vehicles  to  back  out  of 
stalls 

Trip  generations  from  parking  zone  Z 
Parking  arrivals  at  parking  zone  Z 


3.12.3  Reports 

H.3.  PARKING  LOT  TRAVEL  TIMES  AND  DELAYS 


3.13  COORXY 


It  is  necessary  to  know  the  direction  of  flow  on  the  link  (i.e., 
which  end  point  traffic  moves  toward).  This  subroutine  reorganizes  the 
x,y  coordinates  of  any  link,  L,  so  that  X1(L),Y1(L)  represent  the 
upstream  end  of  the  link  and  X2{L),Y2(L)  represent  the  downstream  end 
of  the  link. 


3.14  PLOT! 


3 .14.1  Purpose 

This  subroutine  initializes  values  used  for  plotting  the  volume 
flow  maps.  This  subroutine  also  plots  a  link  or  link  and  zone  map  if 
called  for. 


3.14.2  Data  Description 


Symbol 

EPSLONl 

SCALE 

EPSL0N2 

EPSL0N3 

CHGHT 


XMIN 

XMAX 

YMIN 

YMAX 


Meaning 

Distance  between  volume  lines  in  map  units  (inches) 

The  factor  for  converting  Inches  to  feet  (feet) 

Height  of  hairs  used  to  mark  congestion  (inches) 

Distance  between  hairs  (inches) 

Character  height  of  printer  characters  in  map  units 
(feet) 

X  coordinate  of  lower  left  comer  of  area  to  be  plotted 
X  coordinate  of  upper  right  corner  of  area  to  be  plotted 

Y  coordinate  ol  lower  left  corner  of  area  to  be  plotted 

Y  ((jordltiale  ol  upper  rlp,l>f  corner  ol  area  to  be  plotted 
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3 . 14 . 3  Reports 

A.l.  INPUT  LISTING.  (prints  out  the  XMIN,XMAX,YMIN,YMAX  coordinates) 


3.14.4  Subroutines  Called 

The  following  subroutines  are  called:  PLOT,  PRESCAN,  NUMBER,  ('LOT'S, 
PTRAF,  SYMBOL. 


3,15  PLOTA 

3.15.1  Purpose 

This  subroutine  accumulates  data  that  will  later  be  used  to  make 
summary  plots  of  volumes  on  the  link  network. 


3.15.2  Data  Description 


Symbol 

HRVOL(I) 

DYVOL(I) 

QUEUE(I) 

CAPMAX 

lOVRCAP(l) 


Meaning 

The  sum  for  hourly  volume  plots  for  link  I 
The  sum  for  daily  volume  plots  for  link  I 
The  queue  length  on  link  I 

The  queue  length  that  represents  congestion 

The  overcapacity  links  are  marked  with  a  -1.  for 
a.m.  rush  and  1  for  the  p.m.  rush  hour. 


^3.16  PLOTP 


3.16.1  Purpose 

This  subroutine  causes  a  plot  to  be  written  on  the  output  device, 
writes  a  heading  by  calls  to  Calcomp  routines,  and  calls  PTRAF  to  do  the 
network  plot . 


3.16.2  Data  Description 

Symbol  Meaning 

lUSE(L)  Tells  whether  this  link  L  has  been  plotted 

VOLMAX  The  upper  volume  limit  rounded  off 
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3.16.3  Subroutines  Called 


The  following  subroutines  are  called:  NUMBER,  PLOT,  PTRAF,  SYMBOL. 


3.17  AQAMF 
3.17.1  Purpose 

This  subroutine  accumulates  the  data  results  of  the  BATS  program. 
The  BATS  model  generates  Data  Set  20,  28,  29,  30,  31,  and  32.  The  data 
include  for  each  zone  the  average  speed  in  the  zone  and  the  vehicle 
miles  traveled,  and  for  each  link  the  average  speed  on  the  link  and  the 
number  of  vehicles  driving  the  link.  These  data  must  be  converted  to 
yearly  values  from  the  BATS  results. 


3.17.2  Data  Description 


Symbol 

NUMVA 

NP 

XCENT 

YCENT 

HT 

PL 

AVSPD(Z) 

VMILEM(I) 

SFCS(Z) 

SFHS(Z) 

NHSOAK 

VMILEC(I) 

SFRAMI(Z,1) 

SFRAMI(Z,2) 

NL 

X1(L) 

Y1(L) 

HEIGHT(L) 


Meaning 

Number' of  internal  zones 

Number  of  links  which  surround  a  zone 

Center  of  the  zone — X  coordinate 

Center  of  a  zone — Y  coordinate 

Height 

Length  of  one  side  of  a  zone  (zone  will  be  assumed 
to  be  square  by  AQAM) 

Average  speed  within  the  zone  Z,  miles  per  hour 

Yearly  vehicle  miles  driven  by  vehicles  of  type  I 

Fraction  of  cold  start  vehicles  on  link 

Fraction  of  hot  start  vehicles  on  link 

Number  of  hot  soak  vehicles  parked  in  lot  (recently 
parked  vehicles) 

Yearly  civilian  vehicle  miles  driven  by  vehicles 
of  type  I 

Fraction  of  vehicles  leaving  zone  Z  that  are 
military 

Fraction  of  vehicles  arriving  at  zone  Z  that  are 
military 

Number  of  line  sources  in  network 

X  coordinate  of  line  source  L 

Y  coordinate  of  line  source  L 

Height  above  average  terrain  of  line  source  L 


WIDTH 

Width  of  line  source  (meters) 

X2(L) 

X  coordinate  of  other  end  of  line  source 

L 

Y2(L) 

Y  coordinate  of  other  end  of  line  source 

L 

AVSPD(L) 

Average  speed  on  link  L 

PCCO(L) 

Percent  cold  start  on  link  L 

PCHO(L) 

Percent  hot  start  vehicles  on  link  L 

NCOLDM(K) 

Number  of  cold  military  vehicles  on  link 
type  K 

L  that  are 

NHOTS 

Number  of  hot  soak  vehicles  on  link 

SPRT(Z) 

Sum  of  the  parking  lot  running  times 

SVZ(Z) 

Sum  of  the  parking  lot  arrivals,  zone  Z 

SZV(Z) 

Sum  of  the  parking  lot  departures,  zone  Z 

SCOUNT(I,L) 

Sum  of  the  counts  on  link  L 

3.17.3  Reports 

The  results  of  AQAMF  are  stored  on  Tape  7.  These  are  in  card 
format  and  could  be  printed  or  input  to  the  AQAM  program. 

3. 18  LETTER 

3.18.1  Purpose 

This  subroutine  prints  an  output  page  heading  in  large  block  letters. 

3.18.2  Data  Description 


Symbol 

Meaning 

NYEAR 

The 

date  of  the  run 

LHEAD(l) 

The 

air  base  name,  12  characters 

LHEAD(3) 

The 

run  Identification,  12  characters 

3.18.3  Reports 

Blocked  Run  Title  Page 

3.18.4  Subroutines  Called 

The  following  subroutines  are  called;  TOFC,  CHARAC. 
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■i.l9  TOl'C 
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COUNT ( I, L) 


Count  of  veliicie.s  by  turning  inovoinont ,  Ly()L-,  or 
hot/cold  status  on  link  L 

DIV  The  divisor,  either  1  or  2,  used  to  break  down  trips 

if  they  are  following  more  than  one  route 

3.23  PRESCAN 

3.23.1  Purpose 

This  subroutine  finds  the  maximum  and  minimum  x,y  coordinates  for 
the  air  base. 


3.23.2  Data  Description 

Symbol 
X1(L),X2(L) 

Y1(L),Y2(L) 

3.24  PTRAF 
3.24.1  Purpose 

This  subroutine  plots  the  volumes  associated  with  each  plot  by 
multiple  lines  and  printed  characters.  If  zone  labels  are  called  for  on 
the  plot,  a  call  to  ZLABEL  is  made.  It  also  plots  the  left  boundary 
and  labels  the  corners  with  x,y  coordinates.  Then  it  calls  PPLOT  to 
plot  the  link.  At  the  end  it  plots  the  right  boundary  and  labels  the 
corners . 


Meaning 

X  coordinate  of  link  L 
Y  coordinate  of  link  L 


3.24.2  Data  Description 

Symbol  Meaning 

UPBNDY  Upper  boundary — the  Y  coordinate.  This  limits  the 

size  of  the  plot  when  more  than  one  plot  is  required 
for  a  base. 

LPLTF  Flag  indicating:  0  =  zone  labels  are  not  to  be  Included 

on  any  map;  1  =  zone  labels  are  to  be  included  on  volume 
plots;  2  =  zone  labels  are  to  be  included  on  link  map. 

3.24.3  Reports 

"BOUNDS"  indicates  the  lower  left  X,  upper  right  X,  lower  left  Y 
iiMi  upper  right  Y  bounds  for  the  current  plot. 
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3.24.4  Subroutines  Called 

The  following  subroutines  are  called:  PLOT,  NUMBER,  PPLOT,  SYMBOL, 
Z LABEL. 


3.25  I'PLOT 


3.25. L  Purpose 

This  subroutine  plots  the  links  by  calling  PLOTLK,  starting  with  the 
first  link  and  taking  the  closest  next.  The  lUSE  array  is  changed  from 
"not  used"  =  0,  "used"  =  2  or  "part  used"  =  1.  Subroutine  CLOSEST  is 
used  to  determine  the  next  closest  link. 


3.25.2  Data  Description 

Symbol 
NOTDRW 
IDRAWN 
IPTDRW 
JOKKSCL 

3.25.3  Reports 

If  all  the  links  are  searched  and  no  new  links  are  found,  then 
the  lUSE  array  is  printed  out  without  a  title. 


Meaning 

Link  is  not  drawn  at  all  =  0 
Link  is  completely  drawn  =  2 
Link  is  partly  drawn  =  1 
Link  is  completely  off  scale  =  -1 


3.25.4  Subroutines  Called 


The  following  subroutines  are  called:  PLOTLK,  CLOSEST. 


3.26  CLOSEST 
3.26.1  Purpose 

This  subroutine  finds  the  closest  undrawn  link  and  puts  its  number 
in  NICWLK  to  be  passed  back  to  tiie  calling  program.  It  returns  NEWLK  =  0 
i(  no  closest  link  is  found. 
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3.26.2  Data  Description 


Symbol 
X1(I) 

Y1(I) 

X2(I) 

Y2(I) 

3.26.3  Subroutine  Called 

PDIST  is  the  only  subroutine  called. 

3.27  PDIST 

3.27.1  Purpose 

This  subroutine  finds  the  distance  between  the  current  pen  position 
(PENX.PENY)  and  point  (A,B). 

3.27.2  Data  Description 

Symbol  Meaning 

PENX  The  last  X  position  of  the  pen  in  map  coordinates 

PENY  The  last  Y  position  of  the  pen  in  map  coordinates 

A  X  coordinate  of  point 

B  Y  coordinate  of  point 


Meaning 

Xj^  coordinate  of  link  I 
Yj^  coordinate  of  link  I 
X2  coordinate  of  link  I 
Y^  coordinate  of  link  I 


3.28  PLOTLK 


3. 28. 1  Purpose 

This  subroutine  finds  volume  of  traffic  by  calling  PFNDVOL,  plots 
each  half  street  with  proper  width  lines  by  calling  PLOTLN,  marks  each 
lin*  as  drawn  or  part  drawn,  marks  overcapacity  link*s,  and  labels  links. 


i.28.2  Data  Description 


Symbol 

LK 

LFLC 

TOD 


Meaning 

Link  number 

Flag  Indicating  hourly  (=1)  or  summary  (=2)  plot 
Time  of  day 
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QUEUE (LK) 
CAPMAX 


Length  of  queue  from  which  capacity  measure  made 
Maximum  length  queue  associates  with  congestion 


3.28.3  Subroutines  Called 

The  following  subroutines  are  called;  LABEL,  OVERCAP,  PDISP, 
PFNDVOL,  PLOTLN. 

3.29  PFNDVOL 

3.29.1  Purpose 

This  subroutine  finds  the  volume  of  traffic  on  a  link,  LK,  and 
converts  this  to  a  number  of  lines,  NL.  The  volume  is  saved  in  LABL. 


3.29.2  Data  Description 


Symbo 1 

LK 

NL 

LFLG 

LABL 


Meaning 

Link  number 

Number  of  lines  to  be  plotted  to  show  the  volume 

Flag  for  hourly  or  summary  plot,  LFLG  =  0-no  volume  lines, 
number  links;  LFLG  =  1-hour  volume  lines,  hour  volume 
numbers;  LFLG  =  2-day  volume  lines,  day  volume  numbers 

Label  value  (i.e.,  volume) 


3.30  PDISP 


3.30.1  Purpose 

This  subroutine  displaces  a  plot  line  an  epsilon  distance  to  the 
right  of  the  direction  of  traffic  flow  (which  is  from  XI, Y1  to  X2,Y2) . 


3.30.2  Data  Description 

Symbol 
XI, Y1 
X2,Y2 
EPSLONl 


Meaning 

Upstream  end  of  link 
Downstream  end  of  link 
Distance  between  lines 
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3.31  OVERCAP 


3.31.1  Purpose 

Tills  subroutine  draws  overcapacity  tic  marks  or  lialrs  by  calling 
i'LOTLN  and  RE'iURN  TO  2  wlien  tliere  is  no  line  to  plot  because  it  is  of  I 
the  plot,  i.e.,  beyond  tlie  lower  lei  t  and  upper  riglit  x,y  coord ina les 
specified  for  this  plot. 


3.31.2  Data  Description 


Symbol 

Meaning 

DISTANC 

Distance  XI, Y1  to  X2,Y2 

X1,Y1,X2,Y2 

End  points  of  a  line  to  be  plotted 

A,B 

End  point  of  a  tic  mark 

C,D 

Other  end  of  the  tic  mark — same  as  XI, Y1 

Tl,T2 

End  of  the  tic  mark  corresponding  to  C,D. 
is  offset  for  the  afternoon  period. 

EPSL0N3 

Distance  between  tic  marks 

3  Subroutines 

Called 

The  following  subroutines  are  called:  PDISP2,  PDISP3,  I’hOTLN. 


3.32  LABEL 


3.32.1  Purpose 

This  subroutine  labels  a  line  with  the  volume  count.  The  distance 
between  the  line  and  the  label  is  EPSL0N2. 


3.32.2  Data  Description 


Symbol 

Meaning 

X1,Y1,X2,Y2 

Coordinates  of  the  line  to  be  labeled 

THETA 

Angle  of  the  line 

DIST2 

Distance  to  the  place  where  the  cliaracter  printing 
begins 

NCHAR 

Number  of  characters  to  be  printed 

I  BCD 

Characters  to  be  printed 

ITRKFG 

Flag  telling  if  link  is  completely  plotted 
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HT 

SCALE 


Character  height  in  inches 
Converts  from  feet  to  inches 


3 . 32. 1  Subroutine  Called 

The  following  subroutines  are  called:  CHKBND,  PDISP2,  SYMB02. 

3.33  ZLABEL 

3.33.1  Purpose 

This  subroutine  prints  the  names  or  number  of  the  zones  within  each 

zone. 


3.33.2  Data  Description 
Symbol 


CHGHT 

Character 

height 

XSUM 

Average  of 

the  X 

zone 

YSUM 

Average  of 

the  Y 

zone 

Meaning 

coordinates  of  the  lines  defining  the 
coordinates  of  the  lines  defining  the 


3.33.3  Subroutines  Called 

SYMBOL  is  the  only  subroutine  called. 

3. 34  PLOTLN 

3.34.1  Purpose 

This  subroutine  draws  a  line  between  point  XI, Y1  and  X2,Y2.  If  the 
line  is  off  the  plot  page,  then  it  truncates  the  line  and  sets  a  flag. 

3.34.2  ta  Description 

Symbol  Meaning 

X1,Y1  Coordinates  of  line 

X2,Y2  Coordinates  of  line 


3.34.3  Subroutine  Called 

The  following  subroutines  are  called:  CHKBND,  PDIST,  PLTSCL. 


3.35  PLTSCL 

3.35.1  Purpose 

This  subroutine  scales  the  plot  and  calls  the  CALCOMP  plot  routine 
PLOT  (1969). 


3.35.2  Data  Description 


Symbol 

A 

B 

I 


Meaning 

The  x-coordinate  of  a  point  to  move  the  pen  to 
The  y-coordinate  of  a  point  to  move  the  pen  to 
Indicates  if  the  pen  is  up  or  down  during  the  move 


3.36  CHKBND 
3.36.1  Purpose 

This  subroutine  checks  the  boundaries  of  the  plot  to  see  if  the  line 
extends  beyond  the  boundary,  and  truncates  the  line  if  necessary  by 
calling  XPTMOV  or  YPTMOV. 


3.36.2  Data  Description 


Symbol 

X1,Y1,X2,Y2 

ITRKFG 

DNBNDX 

DNBNDY 

UPBNDX 

UPBNDY 


Meaning 

x,y  coordinates  of  the  line 

Flag  indicating  if  line  is  on  a  plot  window 

Lower  left  X  boundary 

Lower  left  Y  boundary 

Upper  right  X  boundary 

Upper  right  Y  boundary 


3.36.3  Subroutines  Called 

The  following  subroutines  are  called:  XPTMOV,  YPTMOV. 
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3.37 


PTMOVE 


3.37.1  Purpose 

This  subroutine  moves  the  end  points  of  a  line  that  is  out  of  the 
plot  window  along  the  line  to  the  edge  of  the  plot  so  that  the  end 
points  of  the  line  are  within  the  plot  window. 

3.37.2  Data  Uescrlption 

Symbol  Meaning 

XAjYA  Point  which  is  out  of  bounds 

XB,YB  Point  which  may  be  in  bounds 

BOUND  The  X  or  Y  bound 

IB  0  =  test  lower  bound,  1  =  test  upper  bound 

Slope  of  line 
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UTILITY  ROUTINES 


B.l  PROGRAM  NETINT 

The  NETINT  program  (Table  B.l)  uses  the  BATS-LCON  array,  which 
indicates  interconnection  between  links,  to  determine  the  intersections 
in  the  network.  Two  to  four  links  leading  into  each  intersection  are 
specified.  These  link  numbers  are  then  arranged  in  the  order  of  north, 
east,  south,  and  west  approaches.  The  x,y  coordinates  of  each  link  are 
used  to  determine  the  direction  from  which  each  link  approaches  an 
intersection. 

As  input,  the  program  reads  a  header  card  and  a  card  indicating 
the  number  of  links  in  the  network,  followed  by  a  deck  of  Type  2  cards 
in  the  BATS  format.  The  program  prints  an  intersection  number  of  each 
intersection  followed  by  the  link  numbers  of  the  north,  east,  south,  and 
west  approaches.  The  program  also  prints  the  coordinates  and  inter¬ 
section  number  associated  with  each  link. 
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SAMPLE  RUN:  WILLIAMS  AFB 


1.  Introduction 


The  Williams  data  base  has  been  run  through  several  different 
sequences  using  the  BATS  model,  including  runs  for  one  12-hour  time 
period,  for  12  one-hour  time  periods,  for  the  morning  peak  rush  hour 
(15  minute  time  periods  from  7:00  A.M.  to  8:00  A.M.),  for  the  hourly 
period  from  11:00  A.M.  to  12:00  P.M. ,  in  both  descriptive  and  predictive 
modes.  In  this  section  the  11:00  A.M.  to  12:00  P.M.  runs  are  discussed 
and  briefly  analyzed.  A  discussion  of  the  remaining  runs  is  not  within 
the  scope  of  this  document. 


II.  Descriptive  Runs 

The  descriptive  run  is  entitled,  "WILLIAMS  AFB,  11  1300  CAL,  29 
09  79,"  on  the  blocked  heading.  This  heading  is  somewhat  difficult  to 
read  because  the  PPS  (Peripheral  Printing  System)  does  not  overprint 
lines  as  a  line  printer  does. 

The  descriptive  run  simulates  the  VJilllams  AFB  traffic  network 
during  the  peak  noon  hour  from  11:00  A.M.  to  12:00  P.M.  The  vehicle 
volumes  on  each  link  look  quite  reasonable  when  compared  to  field  data 
collected  by  the  research  team.  Thus,  the  run  is  accepted  as  describing 
the  existing  situation  at  the  base  during  a  typical  Friday  noon  rush 
period. 


III.  Predictive  Run 

A  second  run  was  made  to  predict  what  would  happen  if  the  main 
thoroughfares  (D  and  E  Streets)  were  changed  from  two-way  to  one-way 
streets.  A  plot  of  the  link  and  zone  network  is  included  for  reference 
as  Figure  C-1.  This  required  changing  the  input  data  so  that  no  links 
connected  to  any  link  going  west  on  D  Street  or  east  on  E  Street.  Thus, 
the  change  to  the  input  data  was  simple  and  no  changes  were  made  to  the 
intersection  or  zone  definitions. 

The  predictive  run  simulates  the  Williams  AFB  network  during  the 
peak  noon  hour  from  11:00  A.M.  to  12:00  P.M.  The  implementation  of  one¬ 
way  streets  was  tested  to  see  if  improved  traffic  flow  resulted  from 
relieving  congestion  or  by  shortening  turning  movement  time. 

jUi 


IV.  Results 


The  most  general  report  of  network  effectivenss  is  on  the  last  page 
of  the  computer  printout: 

I.l.  NETWORK  SUMMARY  PARAMETERS  FOR  TIME  PERIOD, 

The  results  of  the  descriptive  and  predictive  runs  can  be  summarized 
by  the  following  table: 


Table  C-1 

SUMMARY  OF  PREDICTIVE  AND  DESCRIPTIVE  RUNS 
FOR  WILLIAMS  AFB 


Measured  Factor 

Descriptive 

Predictive 

Total  travel  time  on  network  (vehicle- 
hour) 

71 

75 

Total  running  time  in  parking  zones 
(vehicle-hour) 

73 

73 

Total  vehicle  miles  traveled  on  network 
(vehicle-mile) 

1,603 

1,676 

Total  intersection  delay  (vehicle-hour) 

5 

6 

Total  stops  at  intersection  (vehicle) 

4,818 

5,470 

Total  Intersection  average  queue/lengths 
(meters) 

249 

255 

An  analysis  shows  that  one-way  streets  increase  the  distance 
traveled  by  73  vehicle  miles.  Travel  time  increases  (it  takes  longer 
to  travel  more  miles)  and  intersection  delay  increases  (vehicles  go 
through  more  intersections).  Therefore,  there  is  not  enough  congestion 
during  the  noon  rush  hour  to  offset  the  increase  in  trip  length  caused 
by  one-way  streets.  However,  the  increase  is  not  drastic  and  the  small 
amount  of  congestion  that  was  observed  to  occur  during  the  noon  rush 
hour  could  be  alleviated  by  one-way  streets  and  an  improved  overall 
network  could  result.  A  furtlier  run  of  the  model  could  be  made  by  the 
Williams  AFB  traffic  engineer  to  test  the  hypothesis  for  the  A.M.  and 
P.M.  rush  hours. 
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6.  TRIP  ATTRACTIONS  MODIFIED  BY  GATE  COUNTS  AND  SHIFT 
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D.  TRIP  DISTRIBUTION 

1.  ORIGIN  TO  GATE  AND  GATE  TO  DESTINATION  TRIPS 
( IPFLG(2)>=1 ) 

2.  ORIGIN-DESTINATION  ARRAY  (IPFLO(2):1) 

3.  ORIGIN-DESTINATION  ARRAY  FOR  CIVILIAN  VEHICLE  TRIPS 
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F.  CALIBRATION 

1.  CALIBRATION  FACTORS 
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IIPFLG(3)=2) 

O.  ASSIGNMENT 
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I .  SUMMARY 

1.  NETWORK  SUMMARY  PARAMETERS  FOR  TIME  PERIOD 
POSSIBLE  REPETITION  OF  A  THROUGH  I  FOR  EACH  TIME  PERIOD. 


INTRODUCTION 

THE  U.S.  AIR  FORCE  THROUGH  A  CONTRACTURAL  ARRANGEMENT 
HAS  DEVELOPED  AN  AIR  BASE  MOTOR  VEHICLE  MODEL  THAT  WILL 
SIMULATE  A  BASE  TRAFFIC  NETWORK  USING  AVAILABLE  LAND  USE, 
EMPLOYMENT.  AND  ENGINEERING  DATA.  THE  MODEL  WILL  GRAPHICALLY 
REPRESENT  AIR  BASE  MOTOR  VEHICLE  OPERATION  ON  VOLUME/FLOW 
MAPS,  AND  WILL  OUTPUT  A  FILE  OF  TRAFFIC  FLOWS  FOR  INPUT 
TO  THE  AQAM  (AIR  QUALITY  ASSESSMENT)  MODEL. 1 
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A. 1 .  INPUT  LISTINO:  OF  EACH  DATA  'ARD  -  WITH  MODIFICATIONS  MADE  BY  SUBROUTINE  INPT. 


78/  9/29  WILLIAMS  AFB  11-1300  CALIBRATE 
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5314186.  -0.  6.  11.  -0.  -0.  -0.  -0.  -0.  7.  -0.  -0. 

5324  92.  -0.  5.  1.  -0.  -0.  -0.  -0.  -0.  -0.  -0.  1635. 

5334  50.  -0.  2.  7.  -0.  -0.  2.  -0.  -0.  3.  -0.  -0. 

5341300.  -O.  O.  O.  152.  -0.  -0.  -0.  -0.  -0.  -0.  -0. 

5352592.  -0.  50.  25.  -0.  -0.  -O.  -0.  -0.  -0.  -0.364742. 

5362  9.  -0.  116.  39.  -0.  -0.  -0.  -0.  -0.  3.  -0.  9442. 

5371600.  -0.  -0.  -0.  251.  -0.  -0.  -0.  -0.  -O.  -0.  -O. 

5381200.  -0.  -0.  -0.  70.  -0.  -0.  -0.  -0.  -0.  -0.  -0. 

5392  72.  -0.  12.  14.  -0.  -0.  -0.  -0.  -0.  -0.  -0.  21260. 

5404  70.  -0.  26.  13.  -0.  -0.  -0.  -0.  -0.  -0.  -0.  -O. 

5411300.  -0.  -0.  -0.  102.  -0.  -0.  -0.  -0.  -0.  -O.  -O. 

5422  72.  -0.  26.  29.  -0.  -0.  -0.  -0.  -0.  -0.  -0.  44972. 

5433613.  -0.  26.  204.  -0.  -0.  9.  161.  -0.  4.  -0.  24279. 

5442  50.  -0.  26.  16.  -0.  -0.  -0.  -0.  -O.  -0.  -0.  8512. 

6HOME-WeRK  .50.50*  1  .0005  .0045*  2  .0005  .0045*  3  .0049  .0005*  4  .0049  .0005 
6SHOPP I  NO  *00.103  62 . 48002 . 4800* -0* . 0000* . 0000* -0* . 0000* . 0000* -0« . 0000* . 0000 
6SERVICE  *00.104  7  .3300  .33004  8* . 0000* . 00004  53 . 00006 . 0000* -0* . 0000* . 0000 
6EXTERNAL  .01.015  3  .1565  . 1 648* -0* . 0000* . 0000* -0* . 0000* . 0000* -0* . 0000* . 0000 
6IN0USTRI AL. 50. 502  1 » . 0000* . 0000* -0* . 0000* . 0000* -0* . 0000* . 0000* -0* . 0000* . 0000 
6A0ril NISTER. 50. 506  1  * . 0000* . 0000* -0* . 0000* . 0000* -0* . 0000* . 0000* -0* . 0000* . 0000 
6FLT. LINE  . 50.507  1  * . 0000* . 0000* -0* . 0000* . 0000* -0* . 0000* . 0000* -0* . 0000* . 0000 
6HOME  .50.501  1 • . 0000* . 0000* -0* . 0000* . 0000* -0* . 0000* . 0000* -0* . 0000* . 0000 

6MILITARY  .01.05010  .0003  . 0003* -0* . 0000* . 0000* -0* . 0000* . 0000* -0* . 0000* . OOOO 
7  171  172  371.  366.  -0.  -0.  -0.  -0.  -0.  -0.  -0.  -0. 

9  1.28  1.26  1.28  1.26  1.00  1.00  1.28  1.28  1.28  1.28  1.00  1.00-0.00  600. 

10  14  5171149151229  93231219205207209  44  -0  -0  -0  -0  -0  -0  -0  -0  -0  -0  -0  -0  -0 

10  14  5172150152230  94232220206  42182  -0  -0  -0  -0  -0  -0  -0  -0  -0  -0  -0  -0  -0  -0 
T1*  3.771  T2*  11.931  T2-T1*  8.160 
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INITIALIZATION  UILLIANS  AFB 


t 1-1300  CALIBRATE 

B.1.  ZONE  PARKING  CAPACITIES  AND  TRIP  LENGTHS 


ZONE 

CAPACITY 

LENGTH 

TRAVEL  TINE 

ALTERNATE  TRAVEL  Til 

ns 

-0. 

1603. 

1440. 

1800. 

CH 

-0. 

223. 

760. 

760. 

PH 

-0. 

1604. 

2376. 

2970. 

TM 

-0. 

1603. 

2080. 

2340. 

OL 

-0. 

222. 

648. 

648. 

OT 

-0. 

1604. 

2992. 

2992. 

1 

88. 

227. 

36. 

36. 

2 

lie. 

169. 

26. 

26. 

3 

197. 

243. 

36. 

38. 

4 

170. 

179. 

28. 

28. 

5 

79. 

294. 

40. 

40. 

6 

2%8. 

292. 

46. 

46. 

7 

102. 

298. 

40. 

40. 

a 

30. 

122. 

19. 

19. 

9 

120. 

241  . 

38. 

36. 

10 

97. 

231  . 

36. 

36. 

1 1 

70. 

233. 

36. 

36. 

12 

91  . 

242. 

38. 

38. 

13 

178. 

260. 

41  . 

41  . 

14 

126. 

244. 

38. 

38. 

19 

80. 

232. 

36. 

36. 

16 

38. 

229. 

39. 

39. 

17 

81  . 

137. 

21  . 

21  . 

18 

306. 

268. 

42. 

42. 

19 

20. 

244. 

38. 

38. 

20 

986. 

91  1  . 

80. 

80. 

21 

186. 

248. 

39. 

39. 

22 

168. 

291 . 

39. 

39. 

23 

89. 

140. 

22. 

22. 

24 

199. 

904. 

141  . 

141  . 

29 

82. 

139. 

22. 

22. 

26 

390. 

339. 

92. 

92. 

27 

169. 

290. 

39. 

39. 

26 

64 

239. 

37. 

37. 

23 

217 

397. 

96. 

96. 

30 

29. 

224  . 

39. 

39 

31 

186. 

246. 

39. 

39. 

32 

92. 

903. 

79. 

79. 

33 

90. 

934 

83. 

83. 

34 

300 

590. 

66. 

86 

39 

992. 

462. 

72. 

72. 

36 

9. 

242. 

38. 

38. 

37 

600. 

918. 

81  . 

81  . 

38 

200. 

490 

76. 

76. 

39 

72. 

481  . 

79. 

79. 

40 

70 

649 

101  . 

101  . 

41 

300 

G97. 

103. 

103. 

42 

72. 

142. 

22. 

22. 

43 

613. 

398. 

96. 

96. 

44 

90. 

1923. 

301  . 

301  . 
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•  ••BATS  MODEL  OUTPUT^^*  7B/  9/29  WILLIAMS  AFB  11-1300  CALIBRATE  PERIOD  FROM  HDD.  TO  1200.  HOURS 


<  o  ID  ^  —  N  o  o  ^  ID  N  f>  ID -- r>  •  n  ID  9  o  r>  r»  tA  f>  •  <-0*0 
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0  M  (^  — 


^XOO— -O-OOOOOOOOOOOOOO-OOOOOOMOOOOOO-OOOO-OO-OOOOOOO 

i 

PZlDXZrjK-MO®n«^«90-Mn®«)«N«90-MO^g)®P^090-M0900N®90-MO^ 
C^COC^eO  —-——-—--- -MNMMMMMMMMOOOOOOOOOO^^^^^ 
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TOTAL  ie  4S7  1242  425  CG  119  276  609  426  3G34 


■••BATS  nOOCL  OUTPUT***  70/  9/29  WILLIAMS  APB  11-1300  CALIBRATE  PERIOD  PROM  1100.  TO  1200  HOURS 


^oiDOttno-2P<DNN<D-«>NOMOi-o-tf)5r^oANOOi<v-9)Ain(M(o«-ONM(nr)ir)o»o(MP»cDir) 
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H  ox 

5  S'" 

^  ^3  —  OOOOOOOC0O  —  —  OOOOOO— OOOOOOCMOOOOO  —  OOOOOOOOO  —  oooooo— o 
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ujiiiijujuuu.u.ii.u.k.u. - 0(0-.o<<<co-^oXou:<o"of<o(oc}i2--lS.<oSi^:i 

KOCUA^OO  -•‘--------CMCMCMCMCMCMCMCMCMCMCOnOCOroCOnncOO^iTT^'V 


TCAL  15  ^59  1237  425  111  184  410  607  428  3876 
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—  ir>^'vnoNoa»M<>9iAio(M(0-OM(^Ar>K)AON9><DiD 
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H 
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OOOOOOOOOOOOOOOOOOOOOOOOOOOMOOOOOOOOOOOCWOOnDOQCSlOOO 

^  ^  «  O  —  N 
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—  000000000000000000009)  —  lAOOOeOOOOIOOOOOOOOOOOOOOOOOCO 
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KOOOOOOOOOOOO^-(9NOONOOOOON9»000000000000009»ir>OOCMOOlDC(M 
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o 

z 

zonnionooooooooooooooooooooooooooooooooooooooooooooo 

flC  01  —  N  «  —  O 
UJ  —  — 

I- 

X 

UJ 

UOOOOOOOOOOOOOOOOOCMOOOOO^OOlDOOOOnOOOOOOOOOOOOOOOOOO 
—  NOIOm  N-iMN 

>  O  CM  CM  — 

q: 

UJ 

<0 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


UJZZZZZZ^-I^^^UM«>M«)>>(J>>ZZZ>MM>UJa.^Z>>>UJZ>>>UlM«r>WUI0>U^4.(n 
UJUJUJUJUJUJU.U.U.U.U.Ii.----9}M-«><<<M--«)Z»)U.<»)MMf<Mmv)X  —  —  ZZ  —  MZ-IA.. 

_  W 

DZQZZC^K  —  <Mf)Oir>4Dr»«0lO-NnoO(DNO9IO-CMn^«)<Dr^O^O-<Mn^9)l0NO9>O-Nr>^ 

)-OEOft,)-oo  — —  cycMN<MCM(MCM<vNCMf»Op)nr>r)«or)f)o^nv’r 


•cycMN(MCM(MCM(VNCMf)r)r)nr>r)p)nr>r)o^n^' 


TOTAL  459  >237  425  11 >  184  410  607  428  3876 


AO-A079  555  SRI  INTERNATIONAL  MENLO  PARK  CA  F/6  15/5 

USER  6UI0E  FOR  THE  AIR  FORCE  BASE  AUTOMOTIVE  TRANSPORTATION  SIM— EtC(U) 
SEP  79  R  SANDYS  F08635-76-0-0132 

AFESC/ESL-TR-79-16-V0L-2  ML 


UNCLASSIFIED 


••ATS  MeotL  ts/  •/••  willimo  afs  ii-i90o  calishatc  miso  mm  noo.  Toiaoo. 


I  jbbb0bb,jg|;gjj0-Wbbgbjj--j|j|i3jWjj|j:igr;Kgjb^j8igr;r;;|j^|jb 
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g 

,  i  |f5lgg6'“‘**  —  '‘**S?222PSt5S8588aB«fc8S8s8R*R«fc*»?;5fS3 
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^  W  --  -  --W-  V  ^  ^  « 
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•oo«-ooooooooooooooooo^oooe^ooooooeoooeooooe^ooooooo 
A-»>NON-ooooooo^eeo^oof>-oNOA^o«o^ooo«>eoeoooo0e^oooo*» 

Aoooooooooooooooooooo«>oooo^oooooooooooooooe^oeoooo«» 

W«>ONON«-OOOOOOO«-MO*NO«*«0«>-NYMOM«>-OO-«»-««-ee«-O«Of««-OOO« 

-oooooooooooooooooooooooooeoooooooooeoooeeooooooooe 

oooooooooooooooooooeo»-ooeo^eeoooeeoeooooooo^ooooooo 

•oeoooQooeoeoooooeeooooooeooooooooo-^oeeeoQeoooooeee 
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«««BATS  nOOCL  8UTWT»»»  ft/  •/£•  WILLIAMS  AF»  11-1300  CALIBRATE  PERIOD  FROM  1100.  TOISOO. 

ORI 01  N-OeSTI NATION  ARRAY  (CONTINUED) 
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(•iBATS  nOOCL  OUTPUTax  7S/  9/29  WILLIAMS  AFB  11-1300  CALIBRATE  PERIOD  FROM  1100.  T01200.  HOURS 
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•  ••MTS  MODEL  OOTPUTx*  76/  9/29  WILLIAMS  AFB  11-1300  CALIBRATE  PERIOD  FROM  1100.  Tei200.  HOURS 
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TOTAL  INTERSECTIOW  DELAY  ON  NETWORK  9.(VEH-HRS) 

TOTAL  STO^S  AT  I NTEHSFCT I ONS  4818. (VEH» 

TOTAL  OF  INTERSECTION  AVERAOE  OUEl»F  LENGTHS  ?49.(n) 
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INTRODUCTION 

THE  U.S.  AIR  FORCE  THROUGH  A  CONTRACTURAL  ARRANGEMENT 
HAS  DEVELOPED  AN  AIR  BASE  MOTOR  VEHICLE  MODEL  THAT  WILL 
SIMULATE  A  BASE  TRAFFIC  NETWORK  USING  AVAILABLE  LAND  USE, 
EMPLOYMENT.  AND  ENGINEERING  DATA.  THE  MODEL  WILL  GRAPHICALLY 
REPRESENT  AIR  BASE  MOTOR  VEHICLE  OPERATION  ON  VOLUME/FLOW 
MAPS,  AND  WILL  OUTPUT  A  FILE  OF  TRAFFIC  FLOWS  FOR  INPUT 
TO  THE  AQAM  (AIR  QUALITY  ASSESSMENT)  MODEL . 1 
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AEOYBON.  79/09/16.  SRI  KRONOS/NOS  (0)  WEDNESDAY 


INITIAL  DISTRIBUTION 


AD/DLODL 

1 

HQ  AFISC 

1 

DDC/DDA 

2 

HQ  AUL/LSE  71-249 

1 

HQ  AFSC/DLWM 

1 

HQ  USAFA/Library 

1 

HQ  AFSC/SD 

1 

1  MSEW 

1 

HQ  USAF/LEEV 

1 

OUSDR&E 

1 

HQ  USAF/SGPA 

1 

MTMC/TEA 

1 

OSAF/MIQ 

1 

HQ  AFESC/RDVA 

10 

OSAF/OI 

1 

HQ  AFESC/TST 

2 

AFIT/Library 

1 

SRI  International 

1 

AFIT/DE 

1 

HQ  TAC/DEEV 

1 

NSF 

1 

HQ  SAC/DEPV 

1 

EPA/ORD 

1 

HQ  MAC/DEEE 

1 

USA  Chief,  R&D/EQ 

1 

HQ  ATC/DEPV 

1 

USN  Chief,  R&D/EQ 

1 

HQ  AFLC/DEPV 

1 

OEHL/CC 

1 

USAFRCEAtR 

1 

HQ  AFESC/DEV 

1 

USAFRCE/CR 

1 

USAFSAM/EDE 

1 

USAFRCE/ER 

1 

USAF  Hospital,  Wiesbaden 

1 
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